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HOMOGENIZATION OF INITIAL BOUNDARY VALUE 
PROBLEMS FOR PARABOLIC SYSTEMS 
WITH PERIODIC COEFFICIENTS 

YU. M. MESHKOVA AND T. A. SUSLINA 


Abstract. Let O C R d be a bounded domain of class C 1,1 . In the 
Hilbert space 1/2(0; C"), we consider matrix elliptic second order differ¬ 
ential operators Ad,e and An,s with the Dirichlet or Neumann boundary 
condition on 90, respectively. Here e > 0 is the small parameter. The 
coefficients of the operators are periodic and depend on x/e. The behav¬ 
ior of the operator ] = D,N, for small e is studied. It is shown 

that, for fixed t > 0, the operator e~ A t' et converges in the L 2 -operator 
norm to e~ A i t , as e —¥ 0. Here A® is the effective operator with con¬ 
stant coefficients. For the norm of the difference of the operators e~ A ^ st 
and e _-4 t i a sharp order estimate (of order 0(e)) is obtained. Also, we 
find approximation for the exponential e~ A tM in the (L 2 —»• H 1 )-norm 
with error estimate of order O^s 1 ^ 2 ); in this approximation, a correc¬ 
tor is taken into account. The results are applied to homogenization of 
solutions of initial boundary value problems for parabolic systems. 


Introduction 


The paper concerns homogenization theory of periodic differential opera¬ 
tors (DO’s). A broad literature is devoted to homogenization theory. First 
of all, we mention the books [BaPa] BeLPapl IZhKOj and |Saj . 


0.1. The class of operators. Let O C M d be a bounded domain with the 
boundary of class C 1,1 . We study selfadjoint strongly elliptic second order 
operators Ad,e and An,e acting in the space L2(0; C n ) and depending on the 
small parameter e > 0. Formally, the operators Ad,e and An, £ are given by 
the differential expression 6(D)*g(x/e)6(D) with the Dirichlet or Neumann 
condition on dO, respectively. Here g(x) is a Hermitian (m x m)-matrix- 
valued function assumed to be bounded, positive definite and periodic with 
respect to the lattice T C R d . The operator b( D) is an (m x n)-matrix first 
order DO with constant coefficients. It is assumed that m ^ n. The symbol 
&(£) is subject to some condition ensuring strong ellipticity of the operators 
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under consideration. The simplest example is the scalar elliptic divergence 
type operator — divg(x/<s)V (the acoustics operator). 

Our goal is to find approximations in various operator norms for the 
operator exponential e~ A ^ et , f = D,N, for small e and positive t and to 
apply the results to homogenization of solutions of the parabolic systems. 

0.2. Operator error estimates. Survey. In a series of papers [BSul—3] 
by M. Sh. Birman and T. A. Suslina an operator-theoretic approach to 
homogenization problems was suggested. Operators A e = b(D)*g(x./e)b(D) 
acting in L 2 (M d ;C n ) were studied. In [BSul], it was shown that, as e —>• 0, 
the resolvent (A £ + I) -1 converges in the operator norm in L 2 (M d ;C n ) to 
the resolvent (A 0 +1) -1 of the effective operator A 0 = b(D)*g°b(D). Here 
g° is a constant positive effective matrix. The following error estimate was 
obtained: 

ll(Ar + -0 1 — (A 0 + I) 1 ||L 2 (Kd)_ > .L 2 (Kd) < Ce. (0.1) 

In [BSu3], approximation of the resolvent ( A £ +I )~ 1 in the norm of operators 
acting from L 2 (M d ;C n ) to the Sobolev space H 1 (K d -,C n ) was obtained: 

II (A £ + I) 1 -(A°+I) 1 - eK (e)||L 2 (K d )^-Hi(R d ) ^ Ce. (0.2) 

Here K(s) is the so called corrector; it contains rapidly oscillating factors, 
and so depends on e. We have || sK(£)\\ L2 _y H i = 0(1). 

Estimates (10.11) . (10.21) are order-sharp; the constants in estimates are con¬ 
trolled explicitly in terms of the initial problem data. The results of this 
kind are called operator error estimates. The method of [BSul-3] was based 
on the scaling transformation, the Floquet-Bloch theory and the analytic 
perturbation theory. 

To parabolic problems in W l this method was applied in [Sul-3], [V], 
[VSu], and [M]. In [Sul,2], it was shown that for fixed t > 0 the operator 
exponential e~ Aet converges in the L 2 (M rf ; C n )-operator norm to e~ AOt , as 
e —>• 0; and the following estimate was obtained: 

||e Aet — e A t \\L 2 (S. d )^L 2 (R d ) ^ C £ (t + ff2 ) 1//2 - (°-3) 

In [Su3], approximation of the exponential e~ Ast in the (L 2 —>• )-noriri 
with corrector taken into account was obtained: 

\\e~ Aet - e~ AOt - eJCit-^W^^d^ni^ ^ Ce(t~ l + f~ 1/2 ), t^e 2 . (0.4) 

A different approach to operator error estimates was suggested by 
V. V. Zhikov. In [Zhl,2], [ZhPasl], estimates of the form (10.11) and (10.21) 
were obtained for the acoustics operator and the elasticity operator. In 
[ZhPas2], estimates of the form (10.3 j) and (10.41) were proved for the scalar 
operator A e = — divg(x/e)V. 

Operator error estimates have been also studied for boundary value prob¬ 
lems in a bounded domain; see [ZhPasl], [Grl,2], and [KeLiS]. (A more 
detailed survey can be found in [Su6,8].) 
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For the class of matrix elliptic operators Ad,e and An,e that we consider, 
homogenization problems were studied in recent papers [PSul,2], [Su4-8]. 
In [PSul,2] and [Su4,5], approximations of the operator A~^ e were obtained. 
In [Su6], approximation for the resolvent of An,e was found. The method 
of the cited papers was based on application of the results for the problem 
in M d , introduction of the boundary layer correction term, and a careful 
analysis of this term. Some technical tricks, in particular, using the Steklov 
smoothing, were borrowed from [ZhPasl]. 

For the present paper, the most important are the results of [Su7,8], where 
approximations of the resolvents (Ad,s ~ CO -1 and (An, £ — C-0 -1 with two- 
parametric error estimates (with respect to e and £) were obtained. Let us 
dwell on the main results of [Su7,8]. For ( € C \R+, |£| ^ 1, and 0 < e ^ £o 
(where £q is sufficiently small) we have 

iM t , £ - a) -1 - - ci)~ 1 \\l 2 {o)^l 2 {o) < c^) ( £ |cr 1/2 + £ 2 ), (o.5) 

IIG^t.e-CO 1 ~ (A® - (I) 1 - £Kj(E\ 0 \\l 2 (O)^H 1 (O) 
<C(«( £ >/2| C |-1/4 + e ). <a6) 

Here f = D, N, (j) = arg£; ^4^ is the effective operator given by the expres¬ 
sion 6(D)*g°6(D) with the Dirichlet or Neumann condition; i£f(e;£) is the 
corresponding corrector. These estimates are uniform with respect to 4> in 
any domain of the form {£ € C \ M+, |£| ^ 1, (j)Q ^ (j) ^ 27r — <f> o} with 
arbitrarily small (f>o > 0. Note that for fixed ( estimate (10.51) is of order 
0(e); the order is the same as in estimate (10.11) in M ci . At the same time, 
estimate (10.6 j) is of order O (e 1 / 2 ); the order is worse than in estimate (10.21) 
in R d . This is caused by the boundary influence. 

0.3. Main results. Main results of the paper are approximations of the 
operator exponential f = D,N, in the 7^(0; C n )-operator norm 

and in the norm of operators acting from 7^(0; C n ) to H l (0 ; C n ), with the 
following twoparametric (with respect to e and t) error estimates: 

\\e~ Al tA - e ^||l 2 (0)-kL2(c>) ^ Ce(t + £ 2 )~ l ^ 2 e~ c ^ , t ^ 0, (0.7) 

|| e A.e* — e •'V — e/Cf(t; s)\\l 2 (0)^h 1 (0) ^ Ce 1 ^ 2 t 3 ^ 4 e ; t ^ e 2 , 

(0.8) 

0 < e ^ e 0 , f = D,N. Here is the corresponding corrector. Besides 

(10.81) . we obtain approximation of the operators g(x/e)6(D)e _ “ 4 t,£ t (corre¬ 
sponding to the ,,fluxes”) in the (L 2 —> 7/2)-norm. 

For bounded values of t, the order of estimate (10.71) is sharp and coincides 
with the order of estimate (10.31) for the problem in W l . The order of estimate 
()0.8D deteriorates compared with (10.41) . which is explained by the boundary 
influence. 

Let us explain the origin of the factor e _c t* in (10.71) and (10.81) . In estimates 
(10.31) and (10.41) for the problem in M d such factor is absent, because the point 
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A = 0 is the bottom of the (continuous) spectra of the operators A e and A 0 . 
The spectra of the operators A^ e and are discrete. For f = D these 
operators are positive definite, and one can take cd equal to any positive 
number lying to the left of their spectra. For f = N the number A = 0 is 
an eigenvalue of these operators, but in the difference e~^ N < et — the 

parts of the operators in the kernel Z = Ker An,e = KerAl^ cancel out; 
one can take cjsr equal to any positive number lying to the left of the first 
nonnegative eigenvalues of the operators An, £ and A° N . 

In general case, the corrector contains an auxiliary smoothing 

operator. We distinguish conditions under which it is possible to use the 
standard corrector (without the smoothing operator). 

Moreover, for a strictly interior subdomain O' of the domain O we find 
approximation of the exponential e~ A ^ et in the (L 2 ( 0 ) —» FI 1 ( 0 , ))-operator 
norm with the following error estimate 

||e A i’ st - e -V - e/Cf(f; e)\\ l 2 ( 0 )^h 1 (O') ^ g x 

^ C(S)et~ 1 e~ c ^ t , t ^ e 2 , 0 < e < £q. 

Here the constant depends on 5 = dist {O'; 80}. For bounded values of t 
the order of estimate (10.91) is the same as in estimate (10.41) for the problem 
in M d . This once again shows that deterioration of the order of estimate 
(10.81) is caused by the boundary layer effect. 

0.4. Method. We explain the method of investigation on the example of 
estimate (10.71) in the case of the Dirichlet boundary condition. Let 7 be a 
positively oriented contour in the complex plane enclosing the spectra of the 
operators Ad, £ and A° D . We have 

e ~ A °' et = l e ~ Ct ( AD ’ £ ~ C/) ” 1 dC ’ t > °' 

The exponential of the effective operator satisfies a similar identity. Hence, 
e -A D:e t_ e -A° D t = _A_J e -Ct (G4a£ _ C/) -l _ {A 0 d _ C/) -l } dC (0 . 10) 

Applying the results of [Su 8 ], we obtain approximation of the resolvent for 
( € 7 , and then use representation (10.101) . This leads to (10.71) . Note that the 
character of dependence of the right-hand side of (10.51) on £ for large values 
of |C| is important for us. Approximation with corrector taken into account 
is obtained by the same way. 


0.5. Application of the results to homogenization of solutions of 
initial boundary value problems. The results about approximation of 
the operator exponential e~ A *< st are applied to the study of solutions u e (x, t) 
of the initial boundary value problems for the parabolic equation 


<9u £ (x,f) 


—6(D)*g(x/e)6(D)u e (x, t) + F(x, t), x e O, 0 < t < T, 

( 0 . 11 ) 


dt 





HOMOGENIZATION OF INITIAL BOUNDARY VALUE PROBLEMS 


5 


with the initial condition u E (x, 0) = <^(x) and the Dirichlet or Neu¬ 
mann boundary condition. It is assumed that E L2(0\C n ) and F E 
Lp(( 0, T);L2(0; C n )) with some p. Let uo(x, t) be the solution of the corre¬ 
sponding effective problem. 

We obtain approximations of the solution u e in various norms. In princi¬ 
pal order, we prove estimates of the difference u £ — uo in the L2(0 ; C n )-norm 
for fixed t and in the L p (( 0, T); 1*2(0 ; (C n ))-norm. With corrector taken into 
account, we find approximations of the solution u e in H l (0] C n ) for fixed 
t and in L p ((0, T); H l (0\ C n )). Moreover, for a strictly interior subdomain 
O' we obtain sharper error estimates for approximation of the solution u e 
in H l (0'] C n ) (for fixed t) and in L p ((0, T); iL 1 (0 / ; C n )). 

0.6. Plan of the paper. The paper consists of two chapters and Appen¬ 
dix. In Chapter 1 (§§1-4) we study operators with the Dirichlet boundary 
condition. Chapter 2 (§§5-8) is devoted to operators with the Neumann 
boundary condition. In §1, the class of operators Ad,s is described, the 
effective operator is defined, and the Steklov smoothing operator is intro¬ 
duced. In §2, we formulate the results about the resolvent of the operator 
Ad, e- In §3, the main results of the paper for the case of the Dirichlet 
boundary condition are obtained, namely, approximations for the operator 
in the (L 2 —>• L 2 )- and (L 2 —>• H l (-operator norms are found. Also, 
approximation of the ,,flux” g(x/£)b(D)e~^ D ’ et is obtained. In §4, the re¬ 
sults of §3 about approximation of the operator exponential are applied to 
homogenization of solutions of the first initial boundary value problem for 
the parabolic equation (10.111) . 

The order of presentation in Chapter 2 is similar to that of Chapter 1. 
In §5, the class of operators An,c is described and the effective operator is 
introduced. In §6, the results about the resolvent of the operator An,b are 
given. On the basis of these results, in §7 we obtain approximations for the 
operator exponential e~^ N ’ et and the corresponding fluxes. These are the 
main results of the paper for the case of the Neumann boundary condition. 
In §8, these results are applied to homogenization of solutions of the second 
initial boundary value problem for the equation (10.111) . 

The proofs of the statements about removing the smoothing operator in 
the corrector in the case of additional smoothness of the boundary and in 
the case of estimates in a strictly interior subdomain are given in Appendix 

(§ 9 ). 

0.7. Notation. Let S) and .f)* be complex separable Hilbert spaces. The 
symbols (•, •)/-, and || • ||^ stand for the inner product and the norm in L); the 
symbol || • denotes the norm of a linear continuous operator acting 

from f) to f}*. 

The symbols (•, •) and | • | stand for the inner product and the norm in 
C n ; 1 = l n is the identity (n x n)-matrix. If a is an m x n- matrix, then 
|a| denotes the norm of a viewed as a linear operator from C n to C m . We 
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use the notation x = (x ±,... ,x d ) € M. d , iDj = dj = d/dxj, j = 1 ,d, 
D = — zV = (Di, ... , D d ). 

The L 2 -class of C”-valued functions in a domain O C is denoted by 
L2(0]C n ). The Sobolev classes of C n -valued functions in a domain O C 
R d are denoted by H s (0;C n ). By Hl{P\€. n ) we denote the closure of 
Cq°((D; C n ) in i7 1 (0;C n ). If n = 1, we write simply 1/2(0), H s (0 ), but 
sometimes we use such abbreviated notation also for spaces of vector-valued 
or matrix-valued functions. The symbol L p ((0, T);fj), 1 ^ p ^ oo, stands 
for the Lp-space of ^-valued functions on the interval (0, T). 

We denote M + = [0, oo). By C, C, C, £, , c, and c (possibly, with indices 

and marks) we denote various constants in estimates. 

A brief communication on the results of the present paper is published in 

[MSu]. 

Chapter 1. Homogenization of the first initial boundary value 

PROBLEM FOR PARABOLIC SYSTEMS 

§ 1 . The class of operators Ad, e ■ The effective operator 

1.1. The class of operators. Let T C be a lattice generated by the 
basis ai,..., a<y 

d 

r = {a G : a = rya^, Uj € Z}. 

3 = 1 

Let fl be the elementary cell of the lattice T, i. e., 

d ^ 1 1 

n = {x € R d : X = Tjaj, -- < Tj < -}. 

3 = 1 

We denote |H| = measfb The basis bi,...,b^ in W l dual to the basis 
ai,..., &d is defined by the relations (by a*,) = 5jk- This basis generates a 
lattice T dual to the lattice T. We use the following notation 

ro = — miry |b|, r\ = — diamfh (1.1) 

2 oyber 2 

Below H l (£l) stands for the space of functions in 1L 1 (H) whose T-periodic 
extension to M. d belongs to Hi oc (W l ). For any T-periodic function /(x), 
x S denote 

/ £ ( x ) : = /(e _ix ), ^ > 0. 

Let O C be a bounded domain with the boundary of class C' 1,1 . In 
the space ^(O;C n ), consider the operator Ad, e given formally by the dif¬ 
ferential expression 

A D ,e = 6(D)V(x)6(D) 

with the Dirichlet condition on dO. Here g(x) is a measurable Hermitian 
(m x m)-matrix-valued function (in general, with complex entries) assumed 
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to be periodic with respect to the lattice T, uniformly positive definite and 
bounded. Next, 6(D) is an (m x n)-matrix first order DO of the form 

d 

6(D) = ( 1 - 2 ) 

1=1 

where bi are constant matrices (in general, with complex entries). The 
symbol 6(£) = ^2i=i 6z£z, £ € M. d , corresponds to the operator 6(D). It 
is assumed that m ^ n and 

rank6(£) = n, 0 / £ G R d . (1.3) 

This is equivalent to the existence of the constants «o and a\ such that 

a Q l n ^ b(6)*b(0) ^ a\l n , 9 G § d_1 , 0 < ao ^ < oo. (1.4) 

From (11,21) and (H3D it follows that 

\k\^(h /2 , l = l,...,d. (1.5) 

The precise definition of the operator Ad,e is given in terms of the quadratic 
form 

od |£ [u,u]= f (/(x)6(D)u,6(D)u) dx, u € Hq{0\ C n ). 

Jo 

Under the above assumptions, we have 

Co / |Du| 2 dx ^ aD,e[u, u] ^ c\ / |Du| 2 dx, u G Hq(0; C n ), (1.6) 

Jo ’ Jo 

where Co = ceollff -1 and ci = ckiII^Hloo- It is easy to check these estimates 
extending u by zero to R d \0, using the Fourier transformation, and taking 
O into account. Thus, the form ao i£ is closed in L,2(0;C n ). By the 
Friedrichs inequality and the lower estimate (11.61) . the form ao e is positive 
dehnite: 

aD,e[u,u] ^ c*||u||| 2(CJ) , u € Hl(0;C n ), c* = c o (diam0)~ 2 . (1.7) 

1.2. The effective matrix and its properties. To describe the effective 
operator, we need to introduce the effective matrix g°. Let an (n x m)- 
matrix-valued function A(x) be a weak T-periodic solution of the problem 

6(D)*#(x)(6(D)A(x) + l m ) = 0, ( A(x) dx = 0. (1.8) 

do 

Denote 

ff(x) := ff( x )(6(D)A(x) + l m ). (1.9) 

The effective matrix g° of size mxmis defined by the relation 

9° = l^r 1 / d(x)dx. 

Jn 

It turns out that g° is positive definite. 
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We need the following properties of the effective matrix checked in (BSull 
Chapter 3, Theorem 1.5]. 

Proposition 1.1. The effective matrix satisfies the following estimates 

9 ^9° <9- (1-H) 

Here 

9 = M~ 1 [ g(x)dx, g= f|H | _1 / c/(x) _1 dx 

Jn \ Jn 

If m = n, then the effective matrix g° coincides with g. 

For specific DO’s, estimates (11.111) are known in homogenization theory 
as the Voight-Reuss bracketing. Now we distinguish the cases where one 
of the inequalities in ( 11 . 111 ) becomes an identity (see [ BSull Chapter 3, 
Propositions 1.6 and 1.7]). 

Proposition 1.2. The identity g° =~g is equivalent to the relations 

6 (D)*g fc (x) = 0, k = l,...,m, (1.12) 

where gfc(x), k = 1 ,..., m, are the columns of the matrix gfx). 
Proposition 1.3. The identity g° = g is equivalent to the relations 
lfc(x) = lfc + 6 (D)wfc, 1 £eC™, w k eH\n-C m ), k = 1,... ,m, (1.13) 
where lfc(x), k = 1 ,... ,m, are the columns of the matrix g(x) -1 . 

Obviously, ( 11 . 111 ) implies the following estimates for the norms of the 
matrices g° and (g 0 ) -1 : 

l«°l < Nli„. itoVlsllsTk.. (i.i4) 

1.3. The effective operator. The selfadjoint operator A° D in L 2 ( 0 ;C n ) 
generated by the quadratic form 

a° D [u,u\= f (g°b(D)u,b(D)u) dx, u € Hq (O; C n ), 

Jo 

is called the effective operator for the operator Ad,e From (11.141) it follows 
that the form satisfies the following estimates similar to (11.61) and (11.71) : 

Co / |Du | 2 dx ^ o^[u, u] ^ ci / |Du| 2 dx, u E Hq(0\ C n ), (1.15) 

Jo Jo 

a° D [u,u\^ c4u\\l 2{0 ), u € i7o(0;C n ). (1.16) 

By the condition dO E C 1,1 , the operator A^ is given by the differential 
expression 6 (D)*g° 6 (D) on the domain i7g(0; C n )ni7 2 (0; C n ). The inverse 
operator satisfies 

II (“4-d ) II l 2 {0)^h 2 (o) ^ c • (1-17) 

Here the constant c depends only on ctQ, oq, ||s - 1 ||z,<x,> an d the 

domain O. To justify this fact, we note that the operator &(D)*g° 6 (D) is 
a strongly elliptic matrix operator (with constant coefficients) and refer to 















HOMOGENIZATION OF INITIAL BOUNDARY VALUE PROBLEMS 


9 


the theorems about regularity of solutions of strongly elliptic equations (see, 
e. g., [McL'l. Chapter 4]). 

Remark 1.4. Instead of condition dO G C 1,1 , one can impose the follow¬ 
ing implicit condition on the domain: a bounded domain O with Lipschitz 
boundary is such that estimate (11.171) is satisfied. For such domain the re¬ 
sults of Chapter 1 remain valid (except for the results of Subsection 3.6, 
Theorem 4.3, and Propositions 4.13 and 4.20, where an additional smooth¬ 
ness of the boundary is required). In the case of scalar elliptic operators 
wide conditions on dO sufficient for estimate (11.171) can be found in [KoE] 
and [MaSh, Chapter 7] (in particular, it suffices that dO G C a , a > 3/2). 

1.4. Properties of the matrix-valued function A(x). The matrix¬ 
valued function A(x) belongs to 17 1 (fl). The following estimate was checked 
in [BSu2, (6.28) and Subsection 7.3]: 

IIAIUnn) < M := (|S!|(1 + (2r„)- 2 )ma 0 - 1 || 9 || I ,„|| 9 - 1 |U„) I/2 , (1.18) 

where ro is defined by m- 

We need the following statement proved in [PSu2, Lemma 2.3]. 

Lemma 1.5. Suppose that the matrix-valued, function A(x) is a T-periodic 
solution of the problem (11.81) . Then for any v € Co°(M d ;C m ) we have 

[ |(DA) £ (x)v(x)| 2 dx < Pi [ |v(x)| 2 dx 

J R d J R d 

+ /3 2 e 2 [ |A e (x)| 2 |Dv(x)| 2 dx. 

JR d 

The constants (3 1 and /3 2 depend only on m, d, ao, a\, and ||5 _1 ||l 00 - 

Lemma 1.5 directly implies the following corollary. 

Corollary 1.6. Suppose that the assumptions of Lemma 1.5 are satisfied 
and A G L 00 . Then for any v G H 1 ( R d ; C m ) we have A e v G C n ) and 

ll^ £v ll/I 1 (R d ) ^ 2 II v llz, 2 (R d ) ^JIIAIIioo II v llH 1 (R d ) ‘ 

1.5. The Steklov smoothing. Let S £ be the operator in L 2 (R d ; C m ) given 
by 

(S' e u)(x) = |fl| _1 f u(x — ez)dz (1.19) 

Jn 

and called the Steklov smoothing operator. Note that || ‘S'e || (M d )_>-i^ 2 (iR d ) ^ 1- 

Let s ^ 0. Obviously, D^S/u = S' e D“u for u G LF(M rf ; C m ) and any 
multiindex a such that |ck| ^ s. The operator S £ is a continuous mapping 
of LP(M d ;C m ) to H s ( R d ;C m ), and ||5 £ ||^ (Rd) ^ s(Rd) ^ 1. 

We need the following properties of the operator () 1.19 1) (see [ZhPasll 
Lemmas 1.1 and 1.2] or |PSu2l Propositions 3.1 and 3.2]). 
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Lemma 1.7. Let /(x) be a T-periodic function in such that f £ L 2 (Q). 
Let [f £ ] be the operator of multiplication by the function / £ (x). Then the 
operator [/ £ ]S £ is continuous in L 2 (M. d -,C m ), and 

ll[/ £ ]^ellL 2 (R d HL2(K d ) ^ 1^1 1 ^ 2 |l/lli2(n)- 

Lemma 1.8. Let r\ be defined by (11.11) . Then for any v £ H 1 (K d -, C m ) we 
have 

|](5 £ - /)v|| i 2 (R d } eri||Dv|| L2(Kd) . 


§2. Preliminaries. 

Approximation of the resolvent of the operator Ad,e 

We choose the numbers £o, £i £ (0,1] according to the following condition. 

Condition 2.1. The number £o £ (0,1] is such that the strip 

(dO) £o := {x € M. d : dist {x; 80} < £o} 

can be covered by a finite number of open sets admitting dijfeomorphisms 
of class C 0,1 rectifying the boundary 80. Let £\ = £o(l + ri) -1 , where 
2r\ = diarn f l. 

Obviously, £\ depends only on the domain O and the lattice P. 

Note that Condition 2.1 would be provided only by the condition that 
80 is Lipschitz; we have imposed a stronger restriction 80 £ C 1,1 in order 
to ensure estimate (1LT7P . 


2.1. Approximation of the resolvent of the operator Ad,e i n the 
L 2 (0] C n )-operator norm. We need the results about the behaviour of 
the resolvent ( Ad,e ~ C^) _1 obtained in [Su7,8]. In [Su8 . Theorem 4.1], the 
principal term of approximation of the resolvent was found for ( £ C \ R+ 
such that |C| ^ I- Besides, in [Su8l Theorem 8.1] the (L 2 —> l/ 2 )-estimate for 
the difference of the resolvents ( Ad,e ~ C^) _1 an d ( A° D — (I)~ 1 was obtained 
for f £ C\ [c 0 , oo), where c 0 is a common lower bound of the operators Ad,e 
and A° d . Using (11.71) and (|1.16l) . we take c 0 equal to c* = co(diamO) 2 . 
Below the following notation is systematically used, where <£> £ (0,27r): 


c(ip) 


|sin(£>| 1 , tp £ (0, vr/2) U ( 37 r/ 2 , 27t), 
1, ip £ [vr/2, 37 r/ 2 ]. 


( 2 . 1 ) 


For convenience of further references, the following set of parameters is 
called the initial data: 

m, d , ot o, «i, \\g\\ Laa , \\g Ulooj (2 2 ) 

the parameters of the lattice T; the domain O. 


Theorem 2.2 ( }Su 8 j ). Suppose that O C M. d is a bounded domain with the 
boundary of class C 1,1 . Suppose that the matrix-valued function g(x) and 
DO 6 (D) satisfy the assumptions of Subsection 1 . 1 . Let Ad,e and A Q D be 
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the operators defined in Subsections 1.1, 1.3. Suppose that the number E\ is 
subject to Condition 2.1. 

1°. Let ( € C\M+, C = |C|e^, 4> £ (0, 27r), and |C| ^ 1- Then for 0 < e ^ £i 
we have 

\\{Ad,£ — (I) 1 - (A° D - (I) 1 \\l 2 (o)^l 2 (o) C Ci c (^) 5 (ICI 1//2 e + e 2 ). 

2°. Now, let f £ C \ [c*,oo), where c* = ao\\g~ 1 \\f 1 ^(&\am.O)~ 2 . We put 
C — c* = |C — c* \e 1 ^’, € (0, 2tt), and denote 

(r s = j c W 2 \C - C*\~ 2 , 1C — C*| < 1, 

\c(lp) 2 , |C - c*| ^ 1. 

Then for 0 < e ^ £\ we have 

\\{Ad,£ — (I) 1 - (A° D — (I) 1 \\l 2 (P)^l 2 {0) C C2P*{Cf)e. (2.3) 

The constants C\ and C 2 depend only on the initial data (12.21) . 

2.2. Approximation of the resolvent of the operator Ad, £ in the 
norm of operators acting from L 2 ((D]C n ) to iZ 1 (0;C n ). Approxima¬ 
tion of the resolvent ( Ad,e — C-0 -1 hi the (L 2 —» iZ 1 )-norm was obtained in 
Theorems 4.2 and 8.1 from [Su8| . In order to formulate the result, we intro¬ 
duce the corrector K^(e; C). For this, we fix a linear continuous extension 
operator 


P 0 : H S (Q ; C n ) -A H s (R d - C n ), s ^ 0. (2.4) 


Such a ,,universal” extension operator exists for any bounded domain with 
Lipschitz boundary (see [St] or [R]). Herewith, 


roll 


H s {0)^H s 


c c, 


(s) 
o ’ 


(2.5) 


(s) 

where the constant Cq depends only on s and the domain O. By Ro we 
denote the operator of restriction of functions in onto the domain O. We 
put 

K D (s;C) = Ro [A £ ] S e b(D)Po(A° D - (I)- 1 . (2.6) 

The corrector Kd(e; C) is a continuous mapping of ^(O; C n ) to H l (0\ C n ). 
Indeed, the operator b(D)Po{A° D — C-0” 1 is continuous from L2(0;C n ) to 
iZ 1 (M (i ;C m ). By Lemma 1.7 and relation A € ZZ 1 (H), it is easily seen that 
the operator [A £ ]S £ is continuous from H 1 (M d ;C m ) to Z7 1 (M rf ;C n ). 

Theorem 2.3 ( ]Su8j ). Suppose that the assumptions of Theorem 2.2 are 
satisfied. Let A'£>(e;C) be the operator (12.61) . and letg be the matrix-valued 
function OP . 
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1°. Let C E C\M+, ( = |C|e^ ; f) € (0, 27t), and \(\ ^ 1. Then for 0 < e ^ £i 
we haue 

||(Ad, £ — CO 1 — Md — CO 1 — sL D (e- 0 \\l 2 {O)^h 1 (O) 

< C' 3 c(^) 2 |Cr 1/4 e 1/2 + 0 4 c(</) 4 £, 

h/6(d)(^ D! , - a) -1 - rs e &(D)p 0 (^ - a)- 1 !! L 2 (0)^L 2 {0) 

c c'3c(^) 2 icr 1/ v /2 + c' 4 c(^) 4 e. 

2°. Now, let C E C \ [c*, oo). Then /or 0 < £ ^ £i we have 

II - a)" 1 - - C/)" 1 - etfz>(e;0llMOHHi(O) < ^P*(C)^ 1/2 , 

(2.7) 

||<f &(D)(^, e - Cl)" 1 - rS e b(D)Po(A° D - CO —4 11 l 2 (0)—*l 2 (O) 

C C 5/ o*(C)£ 1/2 . 

The constants C3, C4, C3, C4, C5, and U5 depend only on the initial 
data (12.21) . 

Remark 2.4. In Theorems 2.2 and 2.3, the number c* = 
Q! o||9 ,-1 |Il 1 (diam0)~ 2 can be replaced by any common lower bound c 0 of 
the operators Ad, £ and A° D . Let k > 0 be arbitrarily small number. Due to 
the resolvent convergence, if £ is sufficiently small, one can take c 0 = X® — k, 
where A) 1 is the first eigenvalue of the operator A° D . Under such choice of c 0 
the constants in estimates (E3D, m, and (12.81) will depend on k. 

2.3. The case where A E Loo- It turns out that, under some additional 
assumption on the matrix-valued function A(x), the smoothing operator S £ 
in the corrector can be removed. 

Condition 2.5. Suppose that the T-periodic solution A(x) of problem (11.81) 
is bounded: A E Loo • 

We introduce the operator 

<(e;0 = [A £ ]6(D)(A1 0 d - Cl)- 1 . (2.9) 

By (11.171) . the operator b(D)(A° D — C-0 -1 is continuous from 7/2(0; C n ) to 
H l (0 ; C m ). Applying the extension operator Pq and Corollary 1.6, we check 
that under Condition 2.5 the operator [A £ ] of multiplication by the matrix¬ 
valued function A £ (x) is continuous from H l (0\ C m ) to H l (0\C n ). Hence, 
the operator (12.91) is a continuous mapping of 7/2(0; C n ) to 77 4 (0; (C n ). The 
following result was obtained in [Su8. Theorems 6.1 and 8.3]. 


Theorem 2.6 ( |Su8j ). Suppose that the assumptions of Theorem 2.2 are sat¬ 
isfied. Suppose that the matrix-valued function A(x) satisfies Condition 2.5. 
LetKjj(£-,C) be the operator (12.91) . Letg be the matrix-valued function (11.91) . 
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1°. Let C E C\M+, C = \C\e 1 ^, f € (0, 2ir), and \(\ ^ 1. Then for 0 < e ^ E\ 
we have 

||(Ad, £ — CO 1 - (A° D - Ci) 1 - ^d( £ ; C) 11 z, 2 (c»)—(o) 
C 3 c( 0) 2 |C|- 1/ V /2 + C' 6 c(0) 4 e, 
ii/&(d)m Ac - a) -1 - 5 £ kd)(^d - a)” 1 ii l 2 (o)-l 2( o) 

^ c 3 c( 0 ) 2 icr l/ v / 2 + decide. 

2 °. ./Vow, let C £ C \ [c*, oo). Then for 0 < e ^ ey we have 

II(Ad, £ - a)" 1 - M°d - a)" 1 - eK° D (e- 0h 2 (O)^m(O) < C 7 p.(Oe 1/2 , 

\\g £ b(D)(A D ,e ~ CI)- 1 ~ ¥b(D)(A° D - CI^Wl^l^o) < C 7 p*(C)e l/2 . 

The constants C 3 and C 3 are the same as in Theorem 2.3. The constants 
Cq, Cq, C 7 , and C 7 depend only on the initial data ( 12 . 21 ) and ||AUx,^ - 

In some cases Condition 2.5 is valid automatically. The next statement 
was proved in |BSu3. Lemma 8.7]; see also Remark 8.8 in !>Su3 . 

Proposition 2.7. Condition 2.5 is a fortiori valid if at least one of the 
following assumptions is fulfilled: 

1 ° dimension does not exceed two, i. e., d ^ 2 ; 

2° dimension d is arbitrary, and the operator has the form Ad,e = 
D V(x) D, where the matrix gfx) has real entries ; 

3° dimension d is arbitrary, and the effective matrix g° satisfies g° = g, 
i. e., relations (|1.13j) hold. 

Remark 2.8. Under the assumptions of assertion 2° the norm HAH^ is 
controlled in terms of d, HgH^, ||fl ,_ 1 ||L 00 , and T. Under the assumptions 
of assertion 3° the norm || A|| is controlled in terms of d, m, n, «o, or, 

11 g 11 Loo, \\d 1 11 Loo 5 and the parameters of the lattice T. 

2.4. Estimates in a strictly interior subdomain. In a strictly interior 
subdomain O' of the domain O it is possible to improve error estimates in 
H 1 . The following result was proved in Theorems 7.1 and 8.6 from [Su8| . 

Theorem 2.9 ( [Su8] ). Suppose that the assumptions of Theorem 2.3 are 
satisfied. Let O' be a strictly interior subdomain of the domain O. Denote 
5 :=dist{0';dO}. 

1°. Let ( E C \ M + and |£| ^ 1. Then for 0 < e ^ E\ we have 

\\(Ad,e — (I) 1 - (A° d — CO 1 — £Rd{£‘, C)IIl 2 (0)-)-^ 1 (0') 

c (C 8 r 1 + C' 9 )c(0) 6 e, 

\\fb(-D)(A D , E - CI r 1 ~ g E S e b(B)P o (A 0 D - CI)- 1 \\l 2 (0)^l 2{ o>) 
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2°. Let ( E C\ [c*,oo). Denote p(C) := c(fi)p*(C) + c(-0) 5 / 2 /o*(C) 3 ^ 4 - For 
0 < e ^ £\ we have 

||(Ad, £ - (I) 1 - (A° D - Cl) 1 — eKd{£\ C)llL2(0)-i•^^ 1 (0 , ) 

^ (c w 5~ l p{C) + C'iic(iA) 1/2 P*(0 5/4 ) e, 

\\fb(D)(A D , E - (iy 1 - rS £ b(D)Po(A° D - C0 -1 ||l 2( ohl 2 (O') 

^ (Cior 1 ^) + C n cW 1/2 p*(C) 5/4 ) e- 

The constants Cs, Cg, Cs, Cg, C\g, Cn, C\g, and C\\ depend only on the 
initial data m- 

In the case where A E a similar result holds with a simpler corrector 
Kf)(s',C)- This was shown in Theorems 7.2 and 8.7 from [ ISu8j . 

Theorem 2.10 ( SnHj 'i. Suppose that the assumptions of Theorem 2.9 are 
satisfied. Suppose that the matrix-valued function A(x) satisfies Condition 
2.5. Let K^faC) be the operator (12.91) . 

1°. Let C E C \ M_|_ and |£| ^ 1. Then for 0 < e ^ £i we have 

\\(Ad,e-(I) 1 - (A° d - C-0 1 - £K^ ) (s;C)\\l 2 (0)^h 1 (0') 

^ (C 8 r 1 + Cg)c{fife, 

Ub{p){A D ,e - a)" 1 -g £ b{T>){Al - CI)~ 1 \\l 2 {P)^l 2 ( 0 ') 

^ (C 8 r 1 + Cg)c(4>) 6 £. 

2°. For ( G C \ [c*, oo) and 0 < e ^ E\ we have 

I|(Ad, e —CO 1 - (S\° D — Cl) 1 ~ £K^ ) (£;C)\\l 2 (0)^h 1 (0') 

< (Cio-T^O + C n c(i/0 1/2 p*(C) 5/4 ) e, 

||/6(D)(^ £ - Cl)" 1 -5^(D)(^ - a)- 1 !! L 2 (O)-hL 2 (0') 

^ (CiofT^O + C i nc(V>) 1/2 /9*(0 5/4 ) £• 

The constants Cs, Cs, C\g, and C\g are the same as in Theorem 2.9. The 
constants Cg, Cg, Cu, and C\\ depend only on the initial data (12.2H and 

l|A||x_ 


§3. Approximation of the exponential e A °A 

3.1. The properties of the operator exponential. We start with the 
following simple statement about estimates for the operator exponentials 
e -A D , e t anc i j n various operator norms. 
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Lemma 3.1. Suppose that the assumptions of Theorem 2.2 are satisfied. 


Then for t > 0 and e > 0 we have 

\\e~ AD ’ st \\ L2 (o)^L 2 (o) < e~ Cmt , (3.1) 

lle-^llL 2 (OHtf 1( o) < (co 1 + c7 1 ) 1/2 ^ 1/2 e- c * t/2 , (3.2) 

ll e ADt \\L 2 (0)^L 2 (0) ^ e C **> (3-3) 

ll e ' AjD * 11 z .2 (O)^iY 1 (O) ^ ( c o 1 + c * X ) 1//2 i 1 ^ 2 e c *^ 2 , (3.4) 

ll e ■ A ' d1 |Il 2 (0)^h 2 (0) ^ le c *^ 2 - (3-5) 


Proof. Estimates (j3.1D and <|3.3|> follow directly from (|1.7I) and (11.161) . 
Combining (11.61) and (11.71) . we conclude that 

II u ll at 1 (c>) ^ ( c o 1 + c * 1 )II- / ^z),£ u IIl2(o)’ u e ^o(0 ; < b")- 

Hence, 

|| e -^D, e t|| l^q^hHO) ^ (co 1 + cf 1 ) 1 ^\\A 1 ^ e e~ AD ’ et \\ L ^ 0 )^L 2 (0)- (3-6) 

Next, by (11.71) . 

\\ A D 2 e e ~ AD ’ et \\L 2 (0)^L 2 (0) < sup /Z 1/2 e"^ 

< e-^sup/r^e-^/ 2 < t-i/2 e -^/ 2 . 


Relations (|3.61) and (|3.7f) imply (|3.2D . Similarly, combining the inequality 

UAW^e-^W^o^iO) < t- 1/2 e-“-‘/ 2 , (3.8) 

(11.151) . and (11.161) . we obtain (13.41) . 

By cm 

\\^~ ADt \\L 2 (0)^H 2 (0) ^ c \\-^°D e Ar,t \\L 2 (0)^L 2 (0) ^ c Sup /ie ^ 

^ c * (3.91 

^ ce~ c * t/2 sup /re _Mt/2 < c t _1 e _c * t/2 . 

This implies (13.51) . □ 


3.2. Approximation of the operator exponential in the Z^OjC")- 
operator norm. Our goal in this subsection is to prove the following the¬ 
orem. 

Theorem 3.2. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Then for 0 < e ^ E\ we have 

\\e~ AD ^ - e-^ll l 2 {o)^l 2 (o) < C 12 e(t + £ 2)-l/2 e -c,t/2 ; t ^ 0 (3 . 10 ) 

TTie constant C\ 2 depends only on the initial data (12.21) . 
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Proof. The proof is based on Theorem 2.2 and the representations for expo¬ 
nentials of the operators Ad,e and AP d in terms of the contour integrals of 
the corresponding resolvents. 

We have (see, e. g., [Kal Chapter IX, §1.6]) 

e - A °’*t = -±-Je-«(A D , £ -C I )-i d C, t> 0. (3.11) 

Here 7 is the positively oriented contour in the complex plane enclosing the 
spectrum of Ad,e- By (II.7p and (11.161) . the spectra of the operators Ad,e 
and A^ lie on the ray [c*,oo). Therefore, we take 

7 = {C E C : Im£ ^ 0, ReC = Im ( + c*/ 2 } 
U{CeC:ImC^0, ReC = — Im ( + c*/ 2 }. 


The exponential of the effective operator AP d admits a representation similar 
to (13.111) with the same contour 7 . Hence, 

e --AD,et _ e ~A° D t _ e ~(t ((yl Dj£ _ (■/)”! — (A% — C-0~i) dC. (3.12) 

Using Theorem 2.2, we estimate the difference of the resolvents for £ E 7 
uniformly with respect to arg£. Recall the notation if = arg (£ — c*). Note 
that for £ E 7 and if = n/2 or if = 3n/2 we have |£| = \/5c*/2. We apply 
assertion 2° of Theorem 2.2 for C E 7 such that |£| ^ c := max{l; \/5c*/2}. 
Obviously, if E (it/A, 7ir/A) and p*{C) ^ 2max{l;8c“ 2 } on the contour 7 . 
Therefore, for 0 < e ^ £\ relation (|2.3I) implies that 

I \(A d ,e - CO ' 1 - (A° d - C0 _1 \\l 2 ^l 2 < 2 C 2 max{l; 8 c“ 2 }e 

C C'12 (lCr 1/2 e + £ 2 ) , C€ 7 , |CKc, ^ ^ 

where C' 12 = 2 C *2 maxjl; Sc” 2 }!; 1 / 2 . (The last step is simply coarsening of 
the estimate.) 

For other C E 7 we have | sin <p\ ^ 5 -1 ^ 2 , and, by assertion 1° of Theorem 

2 . 2 , 


ii(Ad,£ - ar 1 - (a° d - ci^Wl^l, c c^ocr 17 ^+oo, 

C € 7, |C| > C, 0 < £ £l, 


(3.14) 


where C" 2 = 5 5 / 2 Ci. As a result, combining (|3. 131) and (13.141) . for 0 < e < E\ 
we have 

iK^-ar'-^-ar'iiL^^^dcr^+e 2 ), (3.15) 

where C '12 = max{C( 2 ; Cf 2 }. 

From (13.121) and (I3~.15l) it follows that 

\\ e ~ AD ^ - e-^H l 2 {0)^l 2 (0) < 27r- 1 Ci2 (V 1/2 r(l/2) + sO" 1 ) e ~ c ^ 2 . 
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(3.16) 


Note that for t ^ e 2 we have e 2 t 1 ^ et 1 ^ 2 . Thus, using the identity 
T(l/ 2 ) = 7r 1 / 2 , we obtain 

We-Ao.ei _ e-^W i2(0) _, i2(0) 47r- 1 / 2 Ci 2 et" 1 / 2 e - c * t/2 

^ Ci 2 £(t + e 2 )~ 1 / 2 e~ c * t/2 , t ^ e 2 , 

where C 12 = 4v / 27r _1 / 2 C'i2. For t ^ e 2 we apply (|3.ip and (13.311 : 

|| £a(0) _> £a(0) < 2e —c *^ < 2^2e(i + e 2 )- 1 /^ - **/ 2 , t < e 2 . 

(3.17) 

Relations (13.161) and (|3.17l) imply the required inequality (13.101) with C 12 = 

max{C’i2; 2\/2}. □ 

3.3. Approximation of the operator e~ AD ’ et in the (L 2 —> 77 1 )-norm. 

We introduce the corrector 

JC D {t-,e) = R 0 [A e ]5 £ 6(D)P 0 e-^ t . (3.18) 

By (13.51) . for t > 0 the operator e~ A n> t is continuous from 7 / 2 ( 0 ;C n ) to 
77 2 (0 ; C n ). Hence, the operator b(~D)Poe~ A °D t is continuous from 7/2(0; C n ) 
to 77 1 (]R ,i ; C m ). As has been already mentioned, the operator [A 2 ]^ is 
continuous from H 1 (R d ;C m ) to 77 1 (M‘ i ;C n ). Thus, the operator /C£>(t;e) is 
a continuous mapping of 7^(0; C n ) to 77 1 (0; C n ). 

Theorem 3.3. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Let /Cr>(t;e) 6 e f/ie operator (13.181) . Let g be the matrix-valued function 
(HD- Then for 0 < e ^ E\ and t > 0 we have 


— eK-d(L, £)\\l 2 (0)^H 1 (0) 

+ et~ 1 )e~ c * t ^ 2 , 

(3.19) 

g £ SebfiD)Poe- A ^\Wo)^ L 2 { o) 

+ st~ 1 )e~ c * t ^ 2 . 

(3.20) 


\\g £ b(D)e~ AD ’ et 
^ 0i 3 (e 1/2 t- 3 ' 

The constants C 13 and C\;$ depend only on the initial data 

Proof. As in the proof of Theorem 3.2, we shall use representations for the 
operators e~ AD ' et and e~ A n t in terms of the integrals of the corresponding 
resolvents over the contour 7 . We have 

e~ AD ’ et - e~ A ° Dt - elC D (t- e) 


1 

2 iri 


7 J e-V {(A d , £ - Cl )- 1 - (A° d - Cl r 1 - eK d {e- 0) dC, 


(3.21) 


where Kjj{e: Q is the operator (12.61) . 

Similarly to (|3. 131) - (13.151) . from Theorem 2.3 it follows that 

II - Cl)-' - (A° d - Cl)-' - cK d (e;C) 

<e 13 (A 2 icr 1/4 


+ £ 


) , C € 7, 0 < e < ei, 


(3.22) 
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with the constant C 13 = max{Cj 3 ; C( 3 }, where C [ 3 = 2 C 5 max{l; 8 c* 2 }c 1 / 4 , 
C{ 3 = max{5C3; 25 C 4 }. Relations (13.211) and (13.221) imply the required 
estimate (13.191) with the constant 6’13 = 27r _1 r(3/4)C'i3. 

In a similar way, inequality (|3.20l) is deduced from the identity 

g £ b(D)e~ AD A - ff e 5 e 6(D)P 0 e-- A ni 

= ~ j (g £ b(D)(A D , e - Cl)- 1 - rS e b(B)Po(A 0 D - CO" 1 ) dQ 

(3.23) 

and Theorem 2.3. □ 

Remark 3.4. Suppose that is the first eigenvalue of the operator A° D . 
Due to the resolvent convergence, for sufficiently small n > 0 there exists a 
number £*, depending on k and on the problem data, such that for 0 < £ ^ £* 
the first eigenvalue Ai j£ of the operator Ad,e satisfies Ai >£ > A^ — k/2. There¬ 
fore, one can transform the contour of integration so that it will intersect 
the real axis at the point A® — k (instead of c*/2). By this way, one can 
obtain estimates of the form (13.101) . (13.191) . and (13.201) with e _c * t//2 replaced 
by e~( a ? _k )*. The constants in such estimates will depend on k. 

3.4. Estimates for small t. Note that for 0 < t < £ 2 it makes no sense 
to apply estimates (13.191) and (13.201) . since it is better to use the following 
simple statement (which is valid, however, for all t> 0 ). 

Proposition 3.5. Suppose that the assumptions of Theorem 2.2 are satis¬ 
fied. Then for t > 0 and £ > 0 we have 

II e~ A oA - e-^H i 2 ( 0 ) ^i (0) < C 1A t- 1/2 e~ cA / 2 , (3.24) 

\\g £ b(fD)e- AD ^\\ L2{0) ^ L2{0 ) < Cut-^e-^ 2 , (3.25) 

|| 5 0 6 (D)e-^ t || L2(o) ^ L2(o) < Cut- 1/2 e~ c * t/2 , (3.26) 

where C\\ = 2a 0 1 ^ 2 || 5 -1 ||^ (l + (diarnO) 2 ) 1 ^ , C 14 = ||</||^. 

Proof. Estimate (|3.24l) follows from (13.21) , (13.41) , and expressions for the con¬ 
stants Co and c*. 

Next, we have 

||/ 6 (D)e _ - 4 £, ' et || l 2 (0)^L 2 (0) < \\9\\L^\\^D^- Al3 ’ et \\L 2 (0)^-L 2 (0)- 

Combining this with (13.71) . we obtain (13.251) . Taking (11.141) and (13.81) into 
account, it is easy to check (13.261) in a similar way. □ 

3.5. The case where A € L^. Now, assume in addition that the matrix¬ 
valued function A(x) satisfies Condition 2.5. In this case it is possible to 
remove the smoothing operator S £ in the corrector and to use the following 
corrector 

JC° D (t-e) = [A £ ]6(D)e-- 4 °n t . 


(3.27) 
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For t > 0 the operator b(D)e~^o t is continuous from 7^(0; C n ) to 
H 1 (0',C m ). As has been already mentioned in Subsection 2.3, under Con¬ 
dition 2.5 the operator [A e ] is continuous from Ff 1 (0;C m ) to H l (0] C n ). 
Hence, the operator Kfi D {t\£) is a continuous mapping of 7^(0; C n ) to 
H'(0:C v ). 

Theorem 3.6. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Let 
Kfi D {t\e) be the operator (13.271) . Then for 0 < e ^ £\ and t > 0 we have 

\\e~ Ar> ’ et - e - - 4 ^ - £tC° D (t] e)\\ L 2 (o)^m( 0 ) < C 15 (e 1 / 2 t~ 3/ 4 + et _1 )e _c * t/2 , 

(3.28) 

\\g £ b(D)e~' AD ’ et - ^b(D)*|U 2(c?) ^z. 2 (O) < Cis^V 3 / 4 + e^V^ 2 , 

(3.29) 

where g is the matrix-valued function CU). The constants C 15 and C 15 
depend only on the initial data ( 12 . 21 ) and ||A||. 

The proof of Theorem 3.6 is similar to that of Theorem 3.3. The difference 
is that instead of Theorem 2.3 one should apply Theorem 2.6. 

3.6. The case of smooth boundary. It is also possible to remove the 
smoothing operator S e in the corrector under the assumption of additional 
smoothness of the boundary. In this subsection we consider the case where 
d ^ 3, because for d ^ 2 Theorem 3.6 is applicable (see Proposition 2.7(1°)). 

Lemma 3.7. Suppose that k ^ 2 is integer. Let O C M. d be a bounded 
domain with the boundary dO of class . Then for t > 0 the operator 

e -A° D t continuous from 7^(0; C n ) to 77 s (0;C n ), 0 < s ^ k, and the 
following estimate is valid: 

\\ e Ar>t \\ l 2 ( 0 )^h s ( 0 ) ^ C s t s / 2 e c * t / 2 ) t > 0. (3.30) 

The constant C s depends only on s, ao, «i, IlffllLoo; II9 ,_1 ||looj an d the do¬ 
main O. 

Proof. It suffices to check estimate (13.301) for integer s € [0, k]: then the 
result for non-integer s follows by interpolation. For s = 0,1,2 estimate 
(I3.30P has been already proved (see Lemma 3.1). 

So, let s be integer and let 2 ^ s ^ k. We apply theorems about regularity 
of solutions of strongly elliptic equations (see, e. g., [McL, Chapter 4]). 
By these theorems, the operator (.d .^,) -1 is continuous from H a (0;C n ) to 
H a+2 (0-C n ) if dO € C a+1,1 . Here er ^ 0 is integer. Note also that the 
operator (^I ^) -1 / 2 is continuous from 7 / 2 ( 0 ; C n ) to 77 x (0; C n ). It follows 
that, under the assumptions of Lemma, for integer s € [2, k\ the operator 
(. A ° d )~ s is continuous from 7^(0; C n ) to 7P(0;C n ). Herewith, 

\\{A° d )~ s ^ 2 \\l 2 (o)^H‘>{o) ^ 


(3.31) 
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where the constant C s depends only on s, ao, cci, HgH^, \\g 1 ||l 00) and the 
domain O. From (J3.31I) it follows that 

\\ e Ar>t \\L 2 ( 0 )-^H s ( 0 ) ^ C s \\(A° D ) s ^ 2 e ADt \\L 2 (0)^L 2 (0) ^ sup x s ! 2 e xt 


x>c r 


^ c q t~ s/2 e~ c * t/2 suva s/2 e~ a/2 ^ C4~ s/2 e ~ c * t/2 , 


sup ck 
a^O 


(3.32) 

□ 


where C s = C s (s/e) s / 2 . 

Using the properties of the matrix-valued function A(x) and the operator 
S e , it is possible to deduce the estimate for the difference of the correc¬ 
tors (13.181) and (|3.27l) from Lemma 3.7 under the assumption of additional 
smoothness of the boundary. 

Lemma 3.8. Let d ^ 3. Let O C be a bounded domain with the boundary 
dO of class C d/1 2,1 if d is even and of class C^ d+1 ^ 2,1 if d is odd. Let K.r){t',e) 
be the operator (13.181) . and let KP D {t\e) be the operator (13.271) . Then 

II^d(L e) — e)||z, 2 (c>)-Kffi(c>) ^ Cd(t 1 +t d / 4 1 ^ 2 )e c * t / 2 ) t > 0. 

(3.33) 

The constant C d depends only on d, m, ao, a\, ||5'||i ;oo , ||g 
domain O. 


-ii 


and the 


(3.34) 

(3.35) 


Lemma 3.8 and Theorem 3.3 imply the following result. 

Theorem 3.9. Suppose that the assumptions of Theorem 2.2 are satisfied, 
and let d ^ 3. Suppose that the domain O satisfies the assumptions of 
Lemma 3.8. LetKP D (t\e) be the operator (13.271) . Let g be the matrix-valued 
function m- Then for 0 < e ^ S\ and t > 0 we have 

\\ e -AD,et _ e -A° D t - £ IC 0 D (t]£)\\ L 2 { 0 ) ^ H i {0) 

^ C * d {e l / 2 t~^ + er 1 + et~ d ^ 2 )e~^ 2 , 

\\ 9 £ b(D)e~ AD ’ et - ?b{B)e- A ^\\ L2(0) ^ L2{ 0) 

< CKe 1/2 t" 3/4 + er 1 + et- d ^ 2 )e~ c ^ 2 . 

The constants and depend only on the initial data m- 

The proofs of Lemma 3.8 and Theorem 3.9 are postponed until Appendix 
(see §9) in order not to overload the main exposition. Clearly, it is convenient 
to apply Theorem 3.9 if t is separated from zero. For small values of t the 
order of the factor (e 1 / 2 ^ -3 / 4 + et~ l + £t~ d ^~ 1 ^ 2 ) grows with dimension. 
This is ,,charge” for removing the smoothing operator. 

Remark 3.10. In Lemma 3.8, instead of the smoothness condition on dO, 
one could impose an implicit condition on the domain: suppose that a 
bounded domain O with Lipschitz boundary is such that estimate (13.311) 
with s = d/2 + 1 is valid. In such domain all the statements of Lemma 3.8 
and Theorem 3.9 are true. 
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3.7. Special cases. The case where the corrector is equal to zero is distin¬ 
guished by the condition g° = g. By Proposition 1.2, in this case relations 
(11.1211 are satisfied, and then the periodic solution A of problem (|1.8I) is 
equal to zero. Applying Theorem 3.3, we obtain the following statement. 

Proposition 3.11. Suppose that the assumptions of Theorem 3.3 are sat¬ 
isfied. Let g° = ~g, i. e. relations (11.121) are satisfied. Then K,£){t\e) = 0, 
and for 0 < e ^ £\ we have 

_ e -A°„, || ij(0) ^„ 1(0) < C 13 (e 1 'V 3 /4 + er')e-^, t > 0. 

If g° = g, by Proposition 1.3, relations (11.131) are satisfied. Then the 
matrix (11,91) is constant: g(x) = g° = g. Besides, by Proposition 2.7(3°), in 
this case Condition 2.5 is satisfied. Applying Theorem 3.6, we arrive at the 
following statement. 

Proposition 3.12. Suppose that the assumptions of Theorem 2.2 are sat¬ 
isfied. Let g° = g, i. e., relations (11.131) are satisfied. Then for 0 < e ^ £\ 
and t > 0 we have 

\\fh(D)e- AD ^ - g°b( ^e-^W^o^o) < C 15 (e 1/2 f 3 / 4 + eT^e^ 2 . 

3.8. Estimates in a strictly interior subdomain. Applying Theo¬ 
rem 2.9, we improve error estimates in a strictly interior subdomain. 

Theorem 3.13. Suppose that the assumptions of Theorem 3.3 are satis¬ 
fied. Let O' be a strictly interior subdomain of the domain O, and let 
6 = dist {O'; dO}. Then for 0 < e ^ e\ and t > 0 we have 

\\e~ AD ' et - e~ A ° Dt - £JC D (t;e)\\ L2 (o)^m(0') < (Cie ^ -1 + C l7 )£t~ l e ~ c * t/2 , 

(3.36) 

\\g £ b(D)e~ AD ’ st - ff £ 5 £ 6(D)P 0 e-^ t || L2(0Hi2(0 , ) 

(C^d - 1 + C 17 )et- 1 e~ cA/2 . 

The constants C\q, C\ 7 , C\ q, and C\ 7 depend only on the initial data (I 2 . 2 |) . 

Proof. The proof is based on application of Theorem 2.9 and identities 
(13.211) . (I3.23p . Besides, we use estimates c(ip) ^ 2 1 / 2 and p*(C) ^ 
2max{l;8cp 2 } for f E 7 , and estimate c(<j)) ^ 5 1 / 2 for f E 7 such that 
|C| > c. These estimates have been checked in the proof of Theorem 3.2. We 
omit the details. □ 

The following result is checked similarly with the help of Theorem 2.10. 

Theorem 3.14. Suppose that the assumptions of Theorem 3.13 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Let 
Kfi D {t;e) be the operator (13.271) . Then for 0 < e ^ £\ and t > 0 we have 

|| e A D , e t _ e A o t — eKfi D (t; £)\\l 2 (P)^hi(0') ^ (CiqS 1 +C\ 7 )et c * 4 / 2 , 

\\g £ b(D)e A °’ £t — <p 6 (D)e ||l 2 (o)^.l 2 (o') ^ (Ci6<5 1 + C'i 7 )et 1 e c *^ 2 . 
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The constants C\§ and C\q are the same as in Theorem 3.13. The constants 
Cyj and C '17 depend only on the initial data m and ||A|| Loo . 

Note that it is possible to remove the smoothing operator S £ from the cor¬ 
rector in estimates of Theorem 3.13 without Condition 2.5. Now we consider 
the case where d 3 (because in the opposite case, by Proposition 2.7(1°), 
Theorem 3.14 is applicable). We know that for t > 0 the operator e~ A ° Dt is 
continuous from L- 2 ( 0 ; C n ) to H 2 (0 ; C n ), and estimate (13.51) holds. Besides, 
the following property of ,,interior regularity” is true: for t > 0 the operator 
e ~A° D t j s con tinuous from L 2 ( 0 ]C n ) to H s (0';C n ) for any integer s ^ 3. 
We have 

\\e~ A ° Dt \\ l 2 ( 0 )^h*( 0 >) < C' s t-nS ~ 2 + t _ 1 ) (s_ 1 ) / 2 e _c * t/2 , t > 0. (3.37) 

The constant C' depends on s , d , «o, an, IMIloo, ||looj and the domain 
O. For scalar parabolic equations the property of ,,interior regularity” was 
obtained in [LaSoUrj (see Chapter 3, § 12). In a similar way it can be checked 
for the operator A ° D ; it is easy to deduce the qualified estimates (|3.37l) using 
that the derivatives D“uo (where uo = with S 7^(0; C n )) sat¬ 

isfy the parabolic equation <9D“uo/3f = — 6 (D)*g° 6 (D)D Q uo. One should 
multiply this equation by x 2 D Q uo and integrate over the cylinder O x ( 0 , t ). 
Here x is a smooth cut-off function which equals zero near the walls and 
the bottom of the cylinder. A standard analysis of the corresponding inte¬ 
gral identity together with the already known estimates of Lemma 3.1 yield 
estimates ()3.37j) . 

The following statement is deduced from (13.371) by using the properties 
of the matrix-valued function A(x) and the operator S £ . 

Lemma 3.15. Suppose that the assumptions of Theorem 3.13 are satisfied, 
and let d ^ 3. Let Kfi D {t] e) be the operator (13.271) . Letr\ be defined by (11.11) . 
Then for 0 < e ^ (4ri) -1 <5 and t > 0 we have 

||/C D (f;e) -K? D (t-,£)\\ L 2 { p)^m{o<) < c ^(t -1 + t~ 1 / 2 { 6~ 2 + t“ 1 ) d/ 4 )e“ c * t/2 . 

(3.38) 

The constant C^ depends only on d, ao, a\, III/ IUooj an d the do¬ 

main O. 

The following result is deduced from Lemma 3.15 and Theorem 3.13. 

Theorem 3.16. Suppose that the assumptions of Theorem 3.13 are satisfied, 
and let d 3. Let Kfi D {t\e) be the operator (13.271) . Then for 0 < e ^ 
rninjei; (4ri) _ 1 <L} and t > 0 we have 

\\e~ AD ’ et _e~ A °ot - £K, Q D {t-,e)\\ L ^ 0 ) ^ Hl { 0 ') < C d eh d {5\ t)e~ c * t/2 , (3.39) 

\\g E b{D)e~ AD ’ st - g £ b{D)e^ A D t \\ l 2 ( 0 )^l 2 ( 0 ') ^ C d eh d { 8 \ t)e~ cA/2 . (3.40) 

Here hd(5;t ) is defined by 

h d (5 ; t) = r 1 ^- 1 + 1) + t “ 1/2 (<T 2 + t~ l ) d ^. 


(3.41) 
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The constants C d and, C d depend only on the initial data Q. 

Proofs of Lemma 3.15 and Theorem 3.16 are postponed until Appendix 
(see §9) in order not to overload the main exposition. Clearly, it is convenient 
to apply Theorem 3.16 if t is separated from zero. For small values of t 
the order of the factor hd(S]t) grows with dimension. This is ,,charge” for 
removing the smoothing operator. 

§4. Homogenization of solutions of the first initial boundary 

VALUE PROBLEM FOR PARABOLIC SYSTEMS 

Now we apply the results of §3 to homogenization of solutions of the first 
initial boundary value problem for parabolic systems. 

4.1. The problem for the homogeneous parabolic equation. We 

start with the homogeneous equation. We are interested in the behaviour 
for small e of the generalized solution u £ of the problem 

= -5(D)*st £ (x)6(D)u e (x,f), x £ 0, t > 0, 

< u £ (x,0) = p(x), xEO, (4.1) 

u £ (-,t)\ d o = o, t > 0, 

where p £ L2(0;C n ). Then u £ (-,t) = e^^Pp. 

Let uo(x, t) be the solution of the so called „homogenized“ problem 

d UO d ( t X,t) = — &(D)*( 7 ° 6 (D)u 0 (x, t), x € O, t> 0, 

< u 0 (x,0) = pipe), xEO, (4.2) 

,u 0 (-,t)|ae> = 0) *>0. 

Then uq(-, t) = e^^^p, and Theorems 3.2 and 3.3 directly imply the fol¬ 
lowing result. 

Theorem 4.1. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Let u £ be the solution of problem m, and let uo be the solution of prob¬ 
lem m, where p £ L2(0 ; C n ). Then for 0 < e ^ £\ and t^ 0 we have 

l|ue(-,t) - u 0 (-,t)|U 2 ( 0 ) ^ C 12 e{t + e 2 ) _1/2 e“ c * t/ - 11^11^(0). 

Suppose that the matrix-valued function A(x) is the V-periodic solution of 
problem USD- Let u 0 (-,f) = PoUo(-,t), where Pq is the extension opera¬ 
tor (12.40 . Let S £ be the smoothing operator (jl.191) . Then for 0 < e ^ £i and 
t > 0 we have 

||u £ (-,f) -u 0 (-,f) - eK £ S £ b(D)u 0 (-,t)\\ H i {0) 

^ C'i 3 (e 1/2 t _3/4 + et“ 1 )e- c * t/2 ||<p|| L2(0) . 

Let g be the matrix-valued function (11.91) . Then the flux p £ := g e b(D)u £ 
admits the following approximation for 0 < e ^ E\ and t > 0: 

Up -g £ S £ b(D)u 0 (-,t)\\ L2(O) ^ Ci 3 (e 1/2 f~ 3/4 + £t~ l )e~ c * t/2 \\p\\ L2{0) . 
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Under the assumption that A € Theorem 3.6 implies the following 
result. 

Theorem 4.2. Suppose that the assumptions of Theorem 4.1 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Then 
for 0 < e ^ £\ and t > 0 we have 

||u £ (-,f) -u 0 (-,f) -eA £ 6(D)u 0 (-,t)||^i (C)) 

< C 15 (£ 1/2 t~ 3/4 + £i _1 ) e~ c * 4/2 11V 7 11 l 2 (O) 5 

IIPeM) -5 £ &(D)u 0 (-,t) \\l 2 (0) < C'i5(e 1/2 i _3/4 + ei" 1 )e _c * t/2 ||^||L 2 (o)- 

In the case of sufficiently smooth boundary from Theorem 3.9 we obtain 
the following result. 

Theorem 4.3. Suppose that the assumptions of Theorem 4.1 are satisfied, 
and let d ^ 3. Suppose that 30 € C d Z 2,1 if d is even and dO € C'( d+1 )/ 2,1 if 
d is odd. Then for 0 < e ^ £\ and t > 0 we have 

||ue(-,t) -u 0 (-,t) -eA £ 6(D)u 0 (-,t)||^i (cl) 

< C^(e 1 / 2 f - 3/4 + EU 1 + £t- d/i - 1 / 2 )e- c * t/ 2 \\v\\ L2(0 ), 

l|Pe(-,*) -F & ( D )uo(-,i)|U 2 (o) 

C^(e 1 / 2 t “ 3/4 + ei _1 + et“ d/ 4 " 1 / 2 )e“ c * t/ 2 ||(p|| i2(0) . 

Special cases are distinguished with the help of Propositions 3.11 and 
3.12. 

Proposition 4.4. Suppose that the assumptions of Theorem 4.1 are satis¬ 
fied. 

1°. If g° = ~g, i. e., relations () 1.1 2 f) are satisfied, then for 0 < £ ^ £\ and 
t > 0 we have 

||u e (-,t) - uo(-,t)|| H i (0) < C 13 (£ 1 / 2 t ~ 3/4 + et“ 1 )e“ c ** /2 ||v?|| L2(0) . 

2°. If g° = g, i. e., relations (11.131) are satisfied, then for 0 < e ^ £\ and 
t > 0 we have 

IIPeM) -PoM)||l 2 ( 0 ) < C'i 5 (e 1/2 f 3/4 + et“ 1 )e" c * t/2 ||<p|| L2(0) , 
where po = g°b( D)uo- 

Approximation of the solutions in a strictly interior subdomain O' of the 
domain O is deduced from Theorem 3.13. 

Theorem 4.5. Suppose that the assumptions of Theorem 4.1 are satisfied. 
Let O' be a strictly interior subdomain of the domain O. Denote 5 := 
dist {O'; dO}. Then for 0 < e ^ £\ and t > 0 we have 

l|ue(-,t) - u 0 (-,i) - £A £ S £ b(D)u 0 (-,t)\\ H i( OI) 

^ (C'le^ -1 + Cn)£t~ 1 e~ c * t ^ 2 \\^p\\ L2 ^ 0 ), 

IIPeM) - 3 e S e &(D)u 0 (-,f)|| L2 ( O /) < {Cif>8~ 1 + Ci 7 )et _1 e _c * t/2 ||<pI|l 2 ( 0 )- 
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In the case where A £ Theorem 3.14 implies the following result. 

Theorem 4.6. Suppose that the assumptions of Theorem 4.5 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Then 
for 0 < e ^ £\ and t > 0 we have 

||u £ (-,f) - u 0 (-,t) - eA £ 6(D)u 0 (-, 011^1(0/) 

^ (Ci6<5 1 + Cn)et l e c * 4// ~ ||V 7 llz. 2 (o)? 

\\Ps(-,t) -g £ b(D)u 0 (-,t)\\ L2 (o') < (Ci6<5 _1 + Ci 7 )et _1 e _c * t / 2 ||v>IU 2 (0)- 

Finally, for d 3 Theorem 3.16 implies the following statement. 

Theorem 4.7. Suppose that the assumptions of Theorem 4.5 are satisfied, 
and let d 3. Then for 0 < e ^ minjei; (4ri)^ 1 d} and t > 0 we have 

||u e (-,i) - u 0 (-,t) - eA £ b(D)u 0 (-,t)\\ H i( O/ ) < C d sh d (6;t)e~ c * t/2 \\(f\\ L ^ 0) , 

\\Pe{-,t) -g £ b(D)u 0 (-,t)\\ L2 ( O ') < C d £h d {5]t)e~ c * tl2 \\^\\ L2{0) . 

Here hfid]t) is given by (13.411) . 

4.2. The problem for nonhomogeneous parabolic equation. Ap¬ 
proximation of solutions in L 2 (0-, C n ). Now, we consider the problem 

dUe gt ’*) = —6(D)*g £ (x)b(D)u e (x, t) + F(x,t), x £ O, 0 < t < T, 

< u e (x, 0) = y?(x), X £ 0, 
u E (-,t)|ao = 0, 0 <t <T, 

(4.3) 

where v? £ L 2 (0;C n ) and F £ Sj p (T) := L p ((0, T); L 2 (0;C n )), 0 < T < oo, 
for some 1 ^ p oo. Then 

u £ (-,t) = + f e -^, s ( t-t) F (.J)dt. (4.4) 

4o 

The corresponding effective problem is given by 

'^£0 =-6(D)V6(D)u 0 (x > t) + F(x > t), x £ 0, 0 < t < T, 

< u 0 (x, 0) = v?(x), X £ 0, 
uo(-,t)|ao = 0, 0 < t < T. 

(4.5) 

Similarly to (14.41) . we have 

u 0 (-,t) = e~ A ° Dt ip(-) + f e~ A ° D ^~^~F(-,t)dt. (4.6) 

Jo 

Subtracting (|4.6I) from (|4.4I) and applying Theorem 3.2, we conclude that 

ll u e(-T) - u o(-T)||l 2 (c>) ^ Ci 2 e(t + e 2 ) 1//2 e cj “ 1111 z. 2 (c?) + C\ 2 £C{£\ t; F), 

(4.7) 

= + \\n-ML 2 { o)dl (4.8) 

Jo 


£(e;t;F) : 
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for 0 < £ ^ £\ and t > 0. 

For 1 < p ^ oo the following result is true. 

Theorem 4.8. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Let u £ be the solution of problem m , and let uo be the solution of problem 
(14.51) for tp G L2(0; C n ) and F € Sj p (T), 0 < T ^ oo, with some 1 < p ^ oo. 
Then for 0 < e ^ E\ and 0 < t < T we have 


||u E (-,t)—u 0 (-,f)|| L2 ( O ) ^ Ci 2 e{t+e 2 ) 1/2 e c * t,2 \\v\\ L 2 (O)+Cp0(e, 
Here 0(e,p) is given by 


6(e,p) = < 


,2-2/p 


1 < p < 2, 


e(| lne| + l) 1 / 2 , p = 2, 
e, 2 < p ^ oo. 

The constant c p depends only on the initial data l\2.2\) and p. 

Proof. By the Holder inequality, from (14.81) we deduce 

£M;F) ^ ||F|| 

r t 

/ o 


(4.9) 


(4.10) 


fie-1) := [\e 2 + r )-p 72 e -c*pV /2 dT _ 

Jo 


(4.11) 


(4.12) 


Here p 1 + (//) 1 = 1 (if p = oo, then p' = 1). For 1 < p < 2, estimating 
the exponential under the integral sign in (|4.12|) by 1, we have 


Lpfis: t) ^ (p'/2 — 1) 1 e 2 p , 1 < p < 2. 


(4.13) 


For p = 2 the integral in (14.121) does not exceed the sum of the integrals 
f ( j (t + e 2 )~ 1 dr and ff° e~ c * T dr. Hence, 

t) ^ 2| In e\ + In 2 + c” 1 ^ max{2; In 2 + c“ 1 }(| lne| + 1). (4.14) 

For 2 < p ^ oo (1 ^ p' < 2) we have 

POO 

Z p fie\t) ^ / r _p, / 2 e _c * p,r / 2 dr = (c*p'/2) p J 2 ~ 1 T(l -p'/2), 2<p^oo. 

Jo 

(4.15) 

Combining (14.71) . (14.111) . and (I4.13p - (l4.15p . we arrive at (14.91) . Herewith, 

'C 12 (p 72-1)-W, 1 < p < 2, 

_ C'i 2 max{ v / 2;(ln2 + cF 1 ) 1 / 2 }, p = 2, 

Cp 1 C , i 2 (c*p , /2) 1/2 -W(r( 1 -p , /2)) 1 /P', 2 < p < oo, 1 ' j 

,C'l2(27r) 1 / 2 C “ 1//2 , p = oo. 


□ 
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4.3. Convergence of the solutions of nonhomogeneous equation 

in L p ((0, T); L2(0; C n )). Now, let 1 ^ p < oo. Using relations (|4.4p and 
(14.61) . positive definiteness of the operators Ad,e and A° D , and the Holder 
inequality, it is easy to check that for F € Sj p (T) the solutions u e and uo 
also belong to Sj p (T). Let us prove the following result. 

Theorem 4.9. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Let u e be the solution of problem (14.31) . and let uo be the solution of problem 
(|4.5I) for ip € L2(0] C n ) and F € Sj p (T), 0 < T ^ oo, with some 1 ^ p < oo. 
Then the solutions u e converge to uo in the $j p (T)-norm, as e —» 0. For 
0 < e ^ £i we have 


K - uo|| flp(T) < c p '®{£,p)\\lp\\ L2{0) + Ci 8 e||F|| %(T) . 
Here ®(e,p) = 9(e,p') is given by 


®(s,p) = 

C 18 = C' 12 (27t) 1 /2 c -V2. Cp/ 
where p _1 + ( p ') _1 = 1. 


< e(| lne| + l) 1 / 2 , 
p 2 tv 

is equal to c p (see 


1 < p < 2, 

P = 2, (4.17) 

2 < p < oo; 

(14. 161) ) with p replaced by p', 


Proof. By (14.71) . 


C{e-t]F) p dt 

(4.18) 

where I p (e; T) = Jq(e 2 + t)~ p / 2 e ~ pc * t / 2 dt (the notation agrees with (14.121) ). 
Using estimates (|4.13l) - (|4.15|) (with p' replaced by p), we see that the first 
summand in the right-hand side of (14.181) does not exceed c p 'Q(E,p)\\<p\\ L2 (Qy 
Now, we consider the second summand in (I4.18p . For p = 1, substituting 
(14.81) and changing the order of integration, we obtain 


|u £ - u 0 


|jjp(T) ^ C 12 £\\<p\\l 2 (0)M £ 'i t ) 1/p + C ^e 



£(e;i;F) dt = / dt ||F(-, t)||z, 2 ( 0 ) / dte c *^ ^/ 2 (e 2 + t - t) ^ 2 
Jo Jt 

< (2/c*) 1/2 r(l/2)||F||* l(r) = (27r/c*) 1/2 ||F|| i j l(T) . 


For 1 < p < oo we apply the Holder inequality to (14.81) : 



\ i/p' 

e -c.(t-t)/2( e 2 + , _ p -l/2 dt \ 



1/p 


(4.19) 
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The integral in the first parenthesis on the right does not exceed (27r/c*) 1//2 . 
Using (14.191) and changing the order of integration, we have 

C(£-t;F) p dt 

< (2vr /c*) p/2p ' £ dt \\F(;t)\\ p L2{<D) £ dte~ c *^/\E 2 + t - i)” 1 / 2 

< (2tt/c*) p/ 2 ||F||^ (t) , 1 < p < oo. 

As a result, the second summand in the right-hand side of (14.181) is estimated 
by CigellFH^^)- ^ 

Remark 4.10. For <p = 0 and F € it is possible to prove also 

convergence of the solutions in o(T). In this case Theorem 4.8 implies 
that 

ll u e — u olUoo(T) ^ Coo£||F|| floo ( T ), 0 < £ ^ £\. 

4.4. The problem for nonhomogeneous parabolic equation. Ap¬ 
proximation of the solutions in H l ((D\ C n ). Now we obtain approxima¬ 
tion of the solution of problem (|4.3|) in C n ) with the help of Theorem 

3.3. Herewith, we have to assume that 2 < p ^ oo. A difficulty arises in ap¬ 
proximation of the integral term in (14.41) . since estimate (13.191) deteriorates 
for small values of t. Assuming that t ^ £ 2 , we represent the integral in (14.41) 
as the sum of two integrals, over the intervals (0, t — £ 2 ) and (t — e 2 , t); the 
first one is estimated with the help of (13.191) , and the second one is estimated 
by using (13.24D . 

We introduce the notation 

w c (-,f) = e~ A D £2 u 0 (-,t - e 2 ), (4.20) 

where uo is the solution of problem (|4.5j) . By (|4.6I) . 

w £ (-,t) = e-^V(0+ / e -^(*-*) F (-,t)dt. (4.21) 

Jo 

Note that w £ (-,t) is the value at the point t of the solution of the first initial 
boundary value problem of the form (14.51) on (0, t) with the new right-hand 
side F e , where F e (-,r) = F(-,r) for 0 < t < t — £ 2 and F e (-,r) = 0 for 
t — £ 2 < t < t. 

Theorem 4.11. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Let u £ be the solution of problem (14.31) . and let uo be the solution of problem 
(14.5j) for (p € L2(0] C n ) and F € Sj p (T), 0 < T ^ oo, with some 2 < p ^ oo. 
Let w £ be defined by (14.201) . Suppose that A(x) is the T-periodic solution 
of problem (11.8[) . Let Pq be the extension operator (12.41) . and let S £ be the 
operator (11.191) . We put w £ = P@w s . Let p £ = g £ &(D)u £j and let g be the 
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(4.22) 

(4.23) 


matrix-valued, function »• Then for 0 < e ^ £\ and e 2 ^ t < T we have 
||u £ (-,t) — u 0 (-,i) — eA £ S e b(D)'w £ (-,t)\\ H i(Q' ) 

< 2Ci 3 £ 1/2 t- 3/4 e- c * t/2 \\ip\\ L2(0) + c p p(e,p) ||F|| flp(t) , 

IIPeM) - 3 £ S £ b(D)we(-,t)\\ L2(0) 

2Ci 3 e 1/2 t- 3/4 e- c * t/2 \\p\\ L2{0) + c p p(E,p)\\F\\sj p{t) . 

Here 

£ 4 ~ 2 / p , 2 < p < 4, 

p{s,p) = < e 1 / 2 (| lne| + l) 3 / 4 , p = 4, (4.24) 

e 1 / 2 , 4 < p ^ oo. 

T/ie constants c p and c p depend only on the initial data m and p. 

Proof. By (11.11) . (14.61) . and (14.211) . we have 

\\u £ (-,t) ~ u 0 (-,t) ~ e^ £ S £ b(D)w £ (-,t)\\ H po) 


+ 

+ 


Ad et _ p 


e 

t-e 2 


- - sKd^e)) ip i 

' H(0) 

-AoAt-t) _ e -A° D (t-t) _ £JCd ^ _ t . j F (. } 


f 

Jo 

Jt (e-ADA*-*) - e -A° D (t-t)^ F(-,t) 


H l {0) 


dt 


H'(0) 


dt. 


(4-25) 

Denote the consecutive summands in the right-hand side of (14.251) by 
J\(t;e;ip), 77 2 (t;£;F), and F). From (|3.19l) and the estimate et~ l ^ 

s 1 / 2 t —3 / 4 , which is valid for t ^ £ 2 , it follows that 


Ji(t]£-ip)^2Ci 3 £ 1/2 t 3/4 e c * t/2 \\<p\\ L2 (0), 0<£^£i. 


(4.26) 


Let us estimate i7 2 (t;e; F). For t E (0,t — £ 2 ) we have t — t ^ e 2 , whence 
e{t — t)~ l ^ £ 1 / 2 (t — f)~ 3 / 4 . Then, by (13.191) . 


J- 2 (t;£;F) ^2C 13 £ 1/2 [ ' (t - t)- 3 / 4 e - c * (t - F) / 2 11F(-,t)|| i2(0) dt, (4.27) 
Jo 

for 0 < £ ^ £i. By the Holder inequality, from (14.271) we obtain 

77 2 (t;£;F) ^ 2C 13 e 1 / 2 ||F|| %(t) I p ,(f, e) 1 ^' , 

H (4.28) 

I p '(t\e) := / r- 3 P'/4 e -c*p'r/2 

(For p = oo and p 7 = 1 the inequality remains true.) 

Now, we estimate I p '{t\e). If 2 < p < 4, then 3p 7 /4 > 1, whence 

POO 

I At ; e) < / t- 3 p'/ 4 dr = j, e 2 ~ 3 P' 12 , 2 < p < 4. (4.29) 

4 e 2 
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If p = 4, then p' = 4/3 and 


' 4/31 


Finally, if 4 < p ^ oo, then 3/7/4 < 1, and 

/•OO 

Ip' (i; e) ^ / t -3p / 4 e _c * p T//2 dr = u p , 4 < p ^ 


Here 


Up = 


(3^/4-I)" 1 , 
max{2; 3(2c*) -1 }, 


3 

d + — ^ i/ 4 (| lnq 
Zc* 

+ !)• 
(4.30) 

4 < p ^ oo. 

(4.31) 

2 < p < 4, 

P = 4, 

(4.32) 

4 < p ^ oo. 

2 < p < oo, 

(4.33) 


As a result, relations (|4.28[) — (14.311) imply that 


where p(s,p) is defined by (14.241) . 

We proceed to the term /73(i;e;F). Taking (j3.24j) into account, we have 

v73(i;e;F) < CiaJ ^t-t)~ 1/2 \\F(-,t)\\ L2{0) dt, 0 < e < 1. (4.34) 

By the Holder inequality, we arrive at 

/73(t;e;F) ^C 14 (l-pV2)- 1/p V- 2 /P||F|| %(i) , 0 < e < 1, (4.35) 

which is true for all 2 < p ^ oo. Note that, if p ^ 4, then e 1 ~ 2 / p ^ e 1 / 2 . 
Combining (14.251) . (14.261) . (14.331) . and (14.351) . we arrive at (14.221) with 


c p = 2C'i 3 Up /p/ + Ci 4 (l — p'/2) 

It remains to check (14.231) . According to (14.41) and (I4.2ip . 

\\Ps(-,t)-g £ S £ b(D)w £ (-,t)\\ L2{0) 


- W /9)- 1 / 'p' 




+ 


+ 


g £ b(D)e~ AD ^ - g e S £ b(B)Poe- A °^ 


<P 


l 2 (o) 


J ^ ( y g e b{D)e~ AD ^ t - t) -g £ Seb^Poe-^-^ F(-,i) 

[ g £ b(D)e- AD ’ s{t - T) F(-,t) 

Jt-e 2 


L 2 (0) 


dt 


L 2 (0) 


(4.36) 

The first two summands in (14.36[) are estimated by using ()3.20l) . and the 
third one is estimated with the help of (13.25[) . (Cf. ()4.26l) — (14.35|) .4 This 

yields (14.231) with the constant c p = 2C\^Vp^ p + Ci 4 (l — p'/2)~ 1 / p '. □ 

Similarly to the proof of Theorem 4.11, under condition A E L <*, Theorem 
3.6 and Proposition 3.5 imply the following result. 
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Theorem 4.12. Suppose that the assumptions of Theorem 4.11 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Then 
for 0 < e ^ £\ and e 2 ^ t < T we have 

l|u £ (-, t) - uo(-,i) -£A e b(D)w £ (-,t)\\ H i io) 

< 2C'i 5 e 1/2 ^ 3 / 4 e _c * t/2 ||^|| L2 ( 0) +c^p(e,p)|| F ||^ (t ), 

\\p £ (-,t) -g £ b(D)w £ {-,t)\\ L2{0) 

< 2C 15 e 1/2 t- 3/4 e- c * t/2 ||^|| L2{0) +c"p( e ,p)||F|| flp(0 . 

T/ie constants c' p and c p depend only on the initial data flZU). ||A|| Loo , and 

p. 

In the case of additional smoothness of the boundary one can apply The¬ 
orem 3.9. However, because of the strong growth of the right-hand sides in 
estimates (I3.34P and (13.351) for small t, we obtain a meaningful result only 
in the threedimensional case and only for p = oo. The proof is similar to 
that of Theorem 4.11. 

Proposition 4.13. Suppose that the assumptions of Theorem 4.11 are sat¬ 
isfied, and let d = 3, p = oo. Suppose that dO € C 2 ’ 1 . Then for 0 < e ^ E\ 
and e 2 ^ t < T we have 

||u £ (-,t) - uoM) -£A e 6(D)w e (-,t)|| H i (0) 

< 2 C 3 ( e V2 r 3/4 + et -5/4^ e- c *t/l\\v\\ L2(0) +c' e l/ 2 || F || floo(t) , 

IIP -g £ b(D)w £ (-,t)\\ L2(0) 

< 2C5 (p/2 t -3/4 + rt -V4) e—/2||^iu s(0) + ? e 1 / 2 ||r|| s „ (t) . 

The constants d and <7 depend only on the initial data (12.21) . 

Now we distinguish special cases. 

Proposition 4.14. Suppose that the assumptions of Theorem 4.11 are sat¬ 
isfied. 

1°. If g° = ~g, i. e., relations (11.121) are satisfied, then for 0 < £ ^ E\ and 
s 2 ^t <T we have 

||u £ (-,f)-uo(-, t)\\ H ^o) < 2C'i3£ 1/2 t" 3/4 e" c * t/2 ||v?|| L2(0) +Cpp(e,p)||F|| flp(t) . 

2°. If g° = g, i. e., relations (I1.13P are satisfied, then for 0 < e ^ £1 and 
£ 2 ^ t < T we have 

||PeM)-PoM)||L 2 ( 0 ) 2C 15 E 1/2 t- 3/4 e- c * t/2 \\<p\\ L2io) +c p p(£,p)\\F\\z p{t ), 

(4.37) 

where po = g°fe(D)uo. The constant d p depends on the initial data (12.2p . 
and also on n and p. 

Proof. Assertion 1° follows directly from Theorem 4.11, since A = 0 if g° = 
V■ 
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Let us prove assertion 2°. By (I4.4(i and (14.611 . 

IIPeM) - Po(-,*)||l 2 (0) 




+ 


9 e b( D)e 




~9°b( D)e 


-4nt 


<P 


l 2 (C) 

5 e 6(D) e -- AD ’ £(t - r) - ff °6(D)e-- 4 o(i-i) 


(4.38) 


FM) 


i 2 (C) 


(it. 


The first term on the right is estimated with the help of Proposition 3.12. 
The integral in the right-hand side of (14.381) can be represented as the sum 
of two integrals, over the intervals (0, t — e 2 ) and (t — £ 2 ,t). The first one 
is estimated by using Proposition 3.12 (cf. (|4.27D - (|4.33p h To estimate the 
second one, we apply (13.2511 and (|3.26l) (cf. (|4.34l) . (14.351) 1. This yields 
estimate (14.37P with d p = 2C\^u p /P +2C' 1 4(1 — p'/2 )~ 1 / p '. □ 

4.5. Approximation of the solutions of nonhomogeneous equation 
in a strictly interior subdomain. Similarly to the proof of Theorem 
4.11, applying Theorem 3.13 and Proposition 3.5, it is easy to obtain the 
following result. 


Theorem 4.15. Suppose that the assumptions of Theorem 4.11 are sat¬ 
isfied. Let O' be a strictly interior subdomain of the domain O. Denote 
5 = dist { O' ; dO}. Then for 0 < £ ^ £\ and e 2 ^t < T we have 

l|ue(-,t) - u 0 (■,*) - eA £ S £ b{D)w £ (-,t)\\ HHO , ) 

< (C'le'J -1 + Cn)£t- l e- c * t/2 \\ip\\ L2 {o) + (M _1 + k' p )a{e,p) ||F|| flp(t) , 

IIP -g e S e b(D)w £ (-,t)\\ L ^ ol) 

< (Ci 6 <5 -1 + Cn)£r 1 e- C * t/2 \\^\\ L2{0 ) + (k p 5~ l + k p )a{£,p)\\F\\^ t y 

Here <r(e,p) is given by 


e 1 2//p , 2 < p < oo, 

e (| lne| + 1), p = oo. 


(4.39) 


The constants k p , k' p , k p , and k' p depend only on the initial data and 

p. 


Finally, in the case where A € L ^ Theorem 3.14 implies the following 
result. 


Theorem 4.16. Suppose that the matrix-valued function A(x) satisfies Con¬ 
dition 2.5. Under the assumptions of Theorem 4.15 for 0 < e ^ £\ and 
e 2 dzt <T we have 

||u e (-,f) -uo(-,t) -eA e &(D)w e (-,i)|| Jf i (0 , ) 

< {Ci & 5~ 1 + Cn)et~ l e~ c * t/2 \\^\\ L2{0) + (k pcT 1 + k' p )a{e,p) ||F|| %W , 
\\Pe(-,t) -g £ b(D)Ws(-fi)\\L 2 (0') 

^ (Ciq6 1 + Cn)et 1 e c * t//2 IIV P llz, 2 (0) + (k p 5 1 + k' p )a(£,p)\\F\\^ p ^ t y 
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The constants k p and k p are the same as in Theorem 4.15. The constants 
k' p and k! p depend only on the initial data m, iiahl., and p. 

Remark 4.17. 1) For p ^ 4 the coefficient at ||F||^ ^ in estimates of The¬ 
orems 4.15 and 4.16 has the order (in e) sharper than the similar coefficient 
in estimates of Theorems 4.11 and 4.12. If 2 < p < 4, these coefficients are 
of the same order so that there is no benefit in transition to a strictly 

interior subdomain. 2) Application of Theorem 3.16 to approximation of the 
solutions of nonhomogeneous equation in the class H l (0'\ C n ) does not lead 
to interesting results, because of the strong growth of the factor hd{S;t ) in 
estimates (I3.39P and (J3.40I) for small t. 


4.6. Approximation of the solutions of nonhomogeneous equation 

in L p ((0,T);fl' 1 (O;C n )). Denote ® P (T) := L p {(0,T)-, H 1 (0;C n )). 


Theorem 4.18. Suppose that the assumptions of Theorem 2.2 are satisfied. 
Let u £ and uo be the solutions of problems (14.3|) and (14.51) . respectively, for 
<p G L2(0; C n ) and F G Sj p (T), 0 < T ^ oo. For e 2 ^ t < T let w £ (-, t) be 
the function defined by (14.201) . For 0 ^ t < e 2 we put w £ (-,t) = 0. Suppose 
that A(x) is theT -periodic solution of problem (11.81) . Let Pq be the extension 
operator m, and let S £ be the operator () 1.19 1) . Denote w £ = Pq w £ . Let 
Pe = g £ b{ D)u e , and let g be the matrix-valued function (USD- 
1°. Let 1 ^ p < 2. Then for 0 < e ^ £\ we have 


l|u e - u 0 - eA e S £ &(D)w £ || 0p(T) ^ K p a(£,p)\\ip\\ L2{0) + Cige^HFH^r), 

(4.40) 

||p £ -5 £ 5 , £ 6(D)w e || ip((0)T);i2(o)) < k p a(£,p)\\<p\\ L2{0) + Ci 9 e 1 / 2 ||F|| %(T) . 
Here a(e,p ) = p(e,p') is defined by 


a(e,p) 


e 1 / 2 , 1 ^ p < 4/3, 

< e 1 / 2 (| Iner| + 1) 3//4 , p = 4/3, 
£ 2 /p- 1 , 4/3 < p <2. 


(4.41) 


The constants C 19 and C 19 depend only on the initial data m- The con¬ 
stants K p and k p depend on the initial data ( 12 . 21 ) and p. 

2°. Let (p = 0 and 1 ^ p ^ 00 . Then for 0 < e ^ E\ we have 

||u £ - u 0 - eA e S’ £ 6 (D)w £ || (gp(T) ^ Cige^F^^), (4.42) 

\\Pe — g 1 ' S e b(D)-W £ \\ Lp ^ 0t T)-L 2 (O)) ^ C , 19^ 1 ^ 2 ||F|| i - :ip ( T ). 

Proof. First, we assume that T ^ e 2 . We have 

||u £ u Q eA £ 5 £ 5(D)w £ || @p(T) ^ Ti(e;p; <p) + 7i(e;p; p>) 

+ T 3 (£-,p-,F) +71(e;p;F) + %(e\p-, F), 


(4.43) 
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where 

Ti(e;p](p) p := 


dt 


(e Ao ’ et - e A ° L,t )‘p 


H 1 


(4.44) 


p 

H 1 


:= [ dt (e A °< et - e ^ - e)C D (t; e)) 

J £ 2 

re 2 rt ^ __ _ 

(e;p;F)P := / dt / - e -^(*-*))F(., t) 

40 Vo 


r 3 

r 4 ( £ ;p;Ff := 


p 

; 

H 1 


P 

; 

H 1 


[ dt f di(e _<A£, - e(t_t) -e~ A o^ - elC D (t-t;e))F(-,t) 

Je 2 Jo 

T 5 (e\p; F)P := [ dt [ £(e~ AD ^-^ - e -^S,(*-*))F(-,i) ? . 

Ve 2 Vt-e 2 Z / 1 

The integral in (14.441) converges for 1 ^ p < 2. Applying (I3.24p . we find 


Ti(e;p-,ip) p ^Cf 4 \\ip\\ p L2(0) j t p / 2 dt = C p lA {\ - p/2) x e 2 p \M P L2(oy 

J 0 

(4.45) 


Next, by (13.191) . 


T 2 {e-p^) p ^{2C 13 ) p e p / 2 \M 

p 


l 2 (o) 


t -3p/4 e -pc*t/2 dL 


(4.46) 


The integral in the right-hand side of (14.461) has been already considered; 
this is I P (T- e) (see (|4.28[) ). We have I p {T ; e) ry for 1 ^ p < 4/3; 
/ 4 / 3 (T; c) ^ ZP 4 (| lne| + 1); / p (T; e) ^ pye 2_3p / 2 for 4/3 < p < 2. Together 
with (14.451) and (14.461) this implies 


Ti(e;p;y)+7^(e;p;y) ^ K p a(e,p)||y|| L2 , 1 < p < 2, (4.47) 

where /%, = Cu(l - p/2)~ 1 / p + 2 Ci 3 iy/ p . 

By the Holder inequality and (|3.24l) . 



For 0 ^ t ^ e 2 we have f^(t — t ) 1 ^ 2 dt ^ 2e. Whence, changing the order 
of integration in the resulting integral, we obtain 


re 2 

r 3 (e;p;F)^Cf 4 (2 £ )^' 

Jo 

sj Cf 4 (2e) 1+p / p '||E 


dt ||F(-, t )||^ 2 

'iip 

N flp(r)- 



dt (t — t) l t 2 


(4.48) 
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The term 7s(e;p; F) is estimated in a similar way: 


%(e ]P ;F)P ^ Cf 4 (2 e ) 1+ ^'||F||^ (T) . (4.49) 

It remains to estimate Ti(e',p;F). By (13.191) and the Holder inequality, 
%(e-p-FY 

^ (2C 13 )P£P/ 2 £ dt^ £ dt(t-t)- 3 / 4 e~ c *^/ 2 \\F(;t)\\L 2 

< (2C 13 ) p £ p / 2 £ dt(J* £ dtit-ty^e-^-V/ 2 \\F(-,t)\\ p L2 

X (/ dt(t- t)- 3 / 4 e - c *h-i)/2 

The last integral in the parenthesis does not exceed (c*/2)^ 1 / 4 r(l/4). 
Changing the order of integration in the resulting integral, we see that 

T 4 (s;p-FY ^ (2 C X3 ) p e p ' 2 ((c*/2) _1/4 r(l/4)) 1+P/P ||F||^ (T) . (4.50) 

Combining ()4.43|) and ()4.47|) - (I4.50I) . we obtain estimate (14.401) (under 
the assumptions of assertion 1°) for T ^ e 2 . Herewith, C\g = 4Ci4 + 
2C'i3(c*/2) _1//4 T(l/4). Obviously, in the case where 0 < T < e 2 estimates 
simplify and rely only on Proposition 3.5. 

Now we assume that ip = 0 and 1 ^ p ^ oo. For 1 ^ p < oo an analog of 
estimate (14.431) holds with 7i(e;p; ip) = 7i(e;p; ip) = 0. Note that estimates 
(14.48D (14.5011 remain true for all 1 ^ p < oo, which implies (14.421) in this 
case. 

For p = oo we have 

||u £ - u 0 - eA £ 5 £ 6(D)w e || (Soo(T) = max{ess sup ||u £ (-,f) - u 0 (-,£)||tfi (o) ; 

tG( 0,e 2 ) 

ess sup ||u £ (•, f) - u 0 (•, t) - £A £ S £ b{F>)w £ (-,t)\\ H i io) }. 
te(e 2 ,T) 

(4.51) 




By Theorem 4.11 for p = oo, 
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where Coo = 2Ci3(c*/2) 1 / 4 r(l/ 4 ) + 2C14. Taking (| 3 . 24 p into account, we 
have 


ess sup ||u £ (-,f) - u 0 (-,t)\\mio) 
te( o,e 2 ) 


= ess sup 
te( o,e 2 ) 


m(o) 


(e~ AD ^ t ~ t) -e- A °o^)F{-,t)dt 

< C 14 ess sup f (t-t) _ 1 / 2 ||F(-,t )|| L2 dt ^ 2eC'i 4 ||F|| i5oo ( T ). 
tefo.e 2 ) do 


(4.53) 


Relations (14.511) (I4.53D imply (|4.42l) for p = 00 . 

The assertion concerning the fluxes p £ can be obtained in a similar way 
with the help of estimates (I3.20p and (|3.25|) . □ 


Under the additional Condition 2.5 on the matrix-valued function A(x) 
the following theorem is true. Its proof is completely analogous to that of 
Theorem 4.18; the difference is that instead of (13.191) and (I3.20p one should 
use (13.281) and (13.291) . 

Theorem 4.19. Suppose that the assumptions of Theorem 4.18 and Con¬ 
dition 2.5 are satisfied. 

1°. Let 1 ^ p < 2. Then for 0 < e ^ £\ we have 
IK - Uo - eA e 6 (D)w £ || (5p(r) ^ x p a(E,p)\\<f\\ L2{0) + C 20 e 1/2 ||F|| flji(r) , 
llPe - if &(D)w £ || Lp (( 0 i T);L 2 (e>)) ^ ^P a ( £ > P) II \\l 2 (O) + ||F|| flp ( T ). 

The constants C 20 and C 20 depend only on the initial data (H 2 D and IIAHi^. 
The constants x p and Tk p depend on the same parameters and p. 

2°. Let if = 0 and 1 ^ p ^ 00 . Then for 0 < e ^ E\ we have 

IK - u 0 - eA £ 6 (D)w £ || 0p(T) < C' 2 oe 1/2 ||F||^ p(T) , 

llPe “ 5 E ^( D )w £ ||lj,(( 0 ,T) ; L 2 (O)) ^ C i 20£ 1 ^ 2 ||F|| flp (T). 

In the case of additional smoothness of the boundary one can apply The¬ 
orem 3.9. As well as in Proposition 4.13, we obtain a meaningful result only 
in the threedimensional case. Its proof is similar to that of Theorem 4.18. 

Proposition 4.20. Suppose that the assumptions of Theorem 4.18 are sat¬ 
isfied, and let d = 3. Suppose that dO E C 2,1 . 

1°. Let 1 ^ p < 4/3. Then for 0 < e ^ we have 

IK - u 0 - £A e 6(D)w e || (gp(r) < x' p s 2/p ~ 3/2 \\<f\\ l 2 (o) + C 2 ie 1/2 ||F|| flij(r) , 
K ~ 5 £& (D) w e||L p ((0,T);L 2 (O)) ^ x' p £ 2/p 3 ^ \\<f\\ L 2 {0) + <?2lK 2 ||F ||fl p (T) • 

The constants C 21 and C 21 depend only on the initial data ( 12 . 21 ) . The con¬ 
stants x! p and Tk' p depend on the same parameters and p. 
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2°. Let ip = 0, and let 1 ^ p ^ oo. Then for 0 < e ^ £i we have 


U £ - Uo - eA e 6(D)w e || Sp(T) < C^i^IlF^^, 
Pe - TKD)we\\l p {(o,t)-,l 2 {0)) ^ C 21 £ 1 / 2 \\F\\ S j p{T) . 


4.7. Approximation of the solutions in L p ((0, T); H l (0'] C n )). By 
analogy with Theorem 4.18, one can check next result with the help of 
Theorem 3.13 and Proposition 3.5. 

Theorem 4.21. Suppose that the assumptions of Theorem 4.18 are sat¬ 
isfied. Let O' be a strictly interior subdomain of the domain O, and let 
5 = dist{0';50}. 

1°. Let 1 ^ p < 2. Then for 0 < e ^ £\ we have 
||u e - u 0 - eA e 5 £ 6(D)w e || ip((0jT); ^i (o , )) 


^ (Aye) 1 + Ap/)r(e,p)||< p|Il 2 (0) + (C' 22 5 1 + C 22 )e (| lne| + l)||F|| fip(T ), 
Pe - g £ S e b( D)w e || £(p((0ir) . £2(O , )) 

^ (Ay <5 1 + A;'/)t(£,p)||<p||l 2 (o) + (C^d 1 + C22MI l ne l + l)ll F llf> P (T)- 


Here r(s,p) = a(e,p') is defined by 



(4.54) 



TTie constants C 22 , C 22 , C 22 , and C 22 depend only on the initial data (|2.2[) . 
T/ie constants Ay, k' pl , Ay, and k' p , are the same as in Theorem 4.15 (with 
p replaced by p'). 

Finally, in the case where A £ Theorem 3.14 and Proposition 3.5 
imply the following result. 

Theorem 4.22. Suppose that the assumptions of Theorem 4.21 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. 

1°. Let 1 ^ p < 2. Then for 0 < e < £i we have 

ll u e — u 0 — eA'6(D)w £ || Z/p (( 0]T ) ;H i( C) /)) 


^ (Ay<5 1 + Ay)r(£,p)||y|| L 2 (e>) + (C 22 S 1 + C 22 )e(\ lne| + l)||F||^ p ( T ), 


Pe - 5 ,£fe (F))w £ || ip ((o,T);L 2 (0')) 

^ (Ay<5 1 + k' p , )t(£,p)\\p\\l 2 ( 0 ) + (C 22 5 1 + C 22 )e(\ lne| + l)||F|| fip ( T ). 
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2°. Let ip = 0, and let 1 ^ p ^ oo. Then for 0 < e ^ £\ we have 
ll u £ — u 0 — eA e 6(D)w £ || Z/p (( 0iT ). H i( C );)) 

^ (C'22^ 1 + C 22 )e(\ l ne l + l)l|F||fl p (T)) 
l|Pe “ #"5 ( D ) w e||Lp((O,T);L 2 (0')) ^ (^22^ 1 + C , 22) £ (l l ne l + 1) ll F lljy(T) ■ 

The constants C 22 and C 22 are the same as in Theorem 4.21. The constants 
C 22 and C 22 depend only on the initial data (12.21) and HAH^. The constants 
kpi, k pl , Ay, and k' p , are the same as in Theorem 4.16 (with p replaced by 

p'). 

Note that application of Theorem 3.16 to approximation of the solutions 
of nonhomogeneous equation in the class L p ((0,T); H 1 (O']C n )) does not 
lead to interesting results because of the strong growth of the factor hd(6] t ) 
in estimates (13.391) and (J3.40I) for small t. 


Chapter 2. Homogenization of the second initial boundary 

VALUE PROBLEM FOR PARABOLIC SYSTEMS 

§5. The class of operators An, s - The effective operator 
5.1. Coercivity. We impose an additional condition on the symbol 6 (£) = 

Ya=i b & for £ G C<i - 

Condition 5.1. The matrix-valued function £>(£), £ E C d , is such that 

rank6(£)=n, 0 / £ E C d . (5.1) 

Note that condition (|5.1D is more restrictive than (11.31) . According to |Ne] 
(see Theorem 7.8 in section 3.7), Condition 5.1 is equivalent to coercivity of 
the form ||i>(D)u||^ 2 ^ in H l (0] C n ). Moreover, the following statement 
holds under the only condition that dO is Lipschitz. 

Proposition 5.2 ( |Ne| ). Condition 5.1 is necessary and sufficient for exis¬ 
tence of constants tfi > 0 and £2 ^ 0 such that the Garding-type inequality 

ll^(D) U ||! 2(0) + £ 2 ||u||| 2(0) >ei||Du||£ a(0) , u eH\0-C n ), (5.2) 

is satisfied. 

Remark 5.3. The constants £1 and £2 depend on the symbol &(£) and 
the domain O, but in the general case it is difficult to control these con¬ 
stants explicitly. However, often for particular operators they can be found. 
Therefore, in what follows we refer to dependence of other constants on 
and £ 2 - 
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5.2. The class of operators. In L 2 ( 0 ;C n ), consider the operator An, £ 
formally given by the differential expression 6 (D)* g £ (x) 6 (D) with the Neu¬ 
mann condition on dO. The precise definition of the operator An,e is given 
in terms of the quadratic form 

ajv, £ [u, u] := f (g e (x) 6 (D)u, 5 (D)u) dx, u € H l (0] C n ). (5.3) 

Jo 

By (11.21) and (11.51) . 

ajv, e [u,u] ^ da;i|| 3 || I/(x) ||Du||| 2(0) , u € H 1 ^; C n ). (5.4) 
From (15.21) it follows that 

F/v,e[u, u] ^ Wg-'Wll (Ci||Du||| 2(0) -£ 2 Ml 2{0) ^ , u eH\0-C n ). 

(5.5) 

Relations (15.31) (15.51) show that the form (15.31) is closed and nonnegative. 

5.3. The effective operator. In L2{0; C n ), consider the operator .4.^ 
generated by the quadratic form 

ajv[u, u] = [ (g°b(D)u, 6 (D)u) dx, u € FI 1 (C ) ; C"). (5.6) 

Jo 

Here g° is the constant effective matrix defined by (jl.lOD . By (|1.14l) , the form 
(|5.6I) satisfies estimates of the form (|5.4j) and (15.51) with the same constants. 

Let i'(x) be the unit external normal vector to dO at the point x E dO. 
Let d [ l be the conormal derivative, i. e., <9°u(x) = b(u(x.))*g°b(V) u(x). 
Since dO E C 1,1 and 6 (D)*g° 6 (D) is a strongly elliptic operator (with con¬ 
stant coefficients), the operator .4.^ can be given by the differential expres¬ 
sion 6 (D)*g° 6 (D) on the domain {u E H 2 (0\ C n ) : 9°u |go = 0}. Herewith, 

\\(A 0 n + I) 1 ||L 2 (O)-kH 2 (0) ^ c °- ( 5 - 7 ) 

The constant c° depends only on the constants and £2 from (15.21) . on crl 
or, Hs'llz.oo? I|s _ 1 ||l cx ,, and ^e domain O. To justify this fact, we refer to 
the results about regularity of the solutions of strongly elliptic systems (see, 
e. g., [ McLl. Chapter 4]). 

The following remark is similar to Remark 1.4. 

Remark 5.4. Instead of condition dO E C 1,1 , one could impose the follow¬ 
ing implicit condition on the domain: a bounded domain O with Lipschitz 
boundary is such that estimate (15.71) is satisfied. For such domain the results 
of Chapter 2 remain valid (except for the results of Subsection 7.6, Theorem 

8.3. and Propositions 8.12 and 8.18, where an additional smoothness of the 
boundary is required). In the case of scalar elliptic operators wide condi¬ 
tions on dO sufficient for estimate (15.71) can be found in [KoE] and [MaSh, 
Chapter 7] (in particular, it suffices that dO E C a , a > 3/2). 
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5.4. Denote 

Z := Ker 6(D) = {z € H 1 (0; <C n ) : 6(D)z = 0}. 

From (15.21) and the compactness of the embedding of i^ 1 (0;C n ) into 
L 2 (C>; C n ) it follows that the dimension of Z is finite. Obviously, Z contains 
the n-dimensional subspace of constant vector-valued functions. Denote 
q = dimZ. We put 

H{0):=L 2 {0- 1 C n )eZ. (5.8) 

Let V be the orthogonal projection of L 2 (0;C n ) onto 71(0). Next, let 
H\{0] C n ) := H\0;C n )nH(0). In other words, 

H±{0] C n ) = {ue H\0 ; C") : (u, z) La(0) = 0, Vz € Z}. 

As shown in [Su 6 l Proposition 9.1], the form ||6(D)u||^ 2 ^ defines a norm in 
the space C n ) equivalent to the standard IL 1 -norm, i. e., there exists 

a constant > 0 such that 

£iH\h H o) < ||6(D)u||| 2(0 ), u e C n ). 

Obviously, 

Ker An,£ = Ker = Z. 

Therefore, the orthogonal decomposition L 2 (0; C n ) = Z 0 71(0) 
the operators An, £ and A° N . By (|5.9D . 

a N,e [ u > u ] > ll9' 1 HZi^ill u ll u € H{(0-X n ), 

4[u,u] ^ ll^llZiCillull^), ueH{(0;C n ). 


(5.9) 

(5.10) 
reduces 

(5.11) 


§6. Preliminaries. 

Approximation of the resolvent of the operator An, £ 


6.1. Approximation of the resolvent of the operator An, £ in the 
L 2 (0; C n )-operator norm. In [Su7,8], approximation of the resolvent of 
the operator (An, £ — C-0 1 was obtained for ( £ C \ [q,,oo), ( A 0- Here 
0 < q, ^ min{/i 2 ,e; H 2 }, where g 2)£ and n® are the first nonzero eigenvalues 
of the operators An, £ and A° N , respectively. (If multiplicities are taken 
into account, they are (q + l)-th eigenvalues.) By (15. lip , one can take 

°b = n^iizici- 

For convenience of further references, the following set of parameters is 
called the initial data: 


m, d, a 0 , ai, || g 1 ||i, 00 ; the constants £i ,£ 2 from (]A 2 J); 

the parameters of the lattice P; the domain O. 


( 6 . 1 ) 


Also, we define the extended data: 

m, n, d, q, a 0 , a u il 9~ i \\L 00 \ 

the constants £], £2 from (15.21) ; the constant £1 from (15.91) ; (6.2) 

the parameters of the lattice F; the domain O. 
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The following result was obtained in Theorems 10.1 and 14.8 of | ISu 8 j . 

Theorem 6.1 ( |Su8j ). Suppose that O C is a bounded domain with the 
boundary of class C 1,1 . Suppose that the matrix-valued function gfx) and DO 
6 (D) satisfy the assumptions of Subsection 1.1. Suppose that Condition 5.1 
is satisfied. Suppose that the number £ 4 € (0,1] is subject to Condition 2 . 1 . 
1°. Let £ € C \ M +; £ = \C,\e 1 ^, and |£| ^ 1. Then for 0 < e =£ £ 4 we have 

\\(AN,,-tur 1 -(A%-or 1 \\L,m^o )« ocw 5 (ici _ 1 / 2 £+p) . (o) 

Here c( f) is defined by (12.11) . The constant C\ depends only on the initial 
data (16.11) , 

2°. Now, let C £ C \ [q,,oo), £ / 0, where q, = ||q _1 ||^Cli. We put 
£ — q, = |£ — c\) |e , d £ (0, 27r) ; and denote 


Pb(C) 


c(??) 2 |£-c b | 2 , |£ — c b | < 1 , 

c(i?) 2 , |£-q,|^l. 


Then for 0 < £ ^ £\ we have 

ll(Av> - C0 _1 - (A% - (iy 1 \\L 2 (0)^L 2 (0) < C2Pb(£)e- (6-4) 

The constant C 2 depends only on the data m - 


6.2. Approximation of the resolvent of the operator An, £ in the 
norm of operators acting from L 2 ( 0 ;C n ) to H l (0', C n ). Approxima¬ 
tion of the operator (An,e — £f )~ 1 with corrector taken into account was 
obtained in ISu81 Theorems 10.2 and 14.81. The corrector is defined similarly 
to ([ 2 TD: 

K N (£-C) = Ro[A e }S £ bfD)Po(A 0 N - (I)- 1 . (6.5) 

The operator Kn(e; £) is a continuous mapping of L 2 ( 0 ; C n ) to H L (0; C n ). 

Theorem 6.2 l [Su 8 ] h Suppose that the assumptions of Theorem 6.1 are 
satisfied. Let A£v(e;£) be the operator (16.51) . where Pq is the extension 
operator m and S £ is the Steklov smoothing operator (11.191) . Let g be the 
matrix-valued function (USD - 

1°. For £ £ C \ |£| ^ 1, and 0 < e =£ £ 4 we have 

II (An,e — (I) 1 - (A° n - £/) 1 — £K n (£- £) ||z. 2 ( 0 )—^iY!(c>) . . 

^ C 3 c(£>) 2 |£r 1 / 4 e 1/2 + C 4 c(£>) 4 e, 

\\g £ b(D)(A N , £ - Cl )' 1 - g £ S £ b(D)P o {A 0 N - a)" 1 1|^(0)^(0) 

^ C 3 c(£>) 2 |£| _ 1 / 4 e 1/2 +C 4 c(</>) 4 £. 


The constants C 3 , C\, C 3 , and C 4 depend only on the initial data (EI]). 

2°. Now, let £ £ C \ [q,,oo), £ / 0. Let V be the orthogonal projection of 
L 2 ( 0 ]C n ) onto the subspace TL(0) defined by (|5.8D . Then for 0 < £ ^ £ 4 












42 


YU. M. MESHKOVA AND T. A. SUSLINA 


we have 

II (* 4 at, e - C-0 -1 - (A% - (I)' 1 - eK N (e\ C,)V\\ L2 (o)^m(0) ^ c 5Pi>(C) £l/2 ; 


|| g £ b(D)(A N , £ - (I r 1 - ?S e b(T>)Po(A° N - (ir^h^^o) 

< C 5 p b (C)e 1/2 . 


( 6 . 8 ) 

(6.9) 


The constants C 5 and C 5 depend only on the data (16.21) . 

Remark 6.3. Instead of ||< 7 ~ J H ^ 1 £ 1 , the number c b in Theorems 6.1 and 
6.2 can be taken equal to any number not exceeding min{// 2 , e ; pfy}- Let 
k > 0 be arbitrarily small number. If e is sufficiently small, we can take 
c \> = P 2 — K ■ Then the constants in estimates will depend on k. 


6.3. The case where A € L^. Suppose that Condition 2.5 is satisfied. 
Then the smoothing operator in the corrector can be removed. We put 

K°n& C) = [A e ]6(D)(^ - Cl)- 1 . (6.10) 

The continuity of the operator iLjy(e;£) from L 2 ( 0 ;C n ) to H 1 (0]C n ) is 
checked similarly to the continuity of the operator (12.91) . 

The following result was proved in Theorems 12.1 and 14.9 of |Su 8 j . 

Theorem 6.4 ( |Su 8 j ). Suppose that the assumptions of Theorem 6.1 are sat¬ 
isfied. Suppose that the matrix-valued function A(x) satisfies Condition 2.5. 
Let (e;C) be the operator (16.101) . Let g be the matrix-valued function 

(HD- 

1°. Let ( € C \ |C| ^ 1. Then for 0 < e ^ E\ we have 

II (>k* - Cl )- 1 - (A% - Clr 1 - eK° N (e-C)\\ L 2 ( 0 )^m( 0 ) 

<: C 3 c(0) 2 |Cr 1/ V /2 +C 6 c(^) 4 e , 

\\fb(D)(A N , £ - Cl )- 1 - g £ b(n)(A - a) -1 II l 2 ( 0 )^l 2 ( 0 ) 

< C3c(^) 2 |cr 1/4 e 1/2 + C 6 c((j)) 4 e. 

The constants C 3 and C 3 are the same as in Theorem 6.2. The constants Cq 
and Cq depend only on the initial data (ED and ||Allioo • 

2°. Now, let C E C \ [c b , 00 ), C A 0. Then for 0 < e ^ £\ we have 

II {An,e — Cl) 1 — {A% — Cl) 1 — £ Kn( £ i C)IIl 2 (0)->H' 1 (0) ^ C7p b (C)£^ 2 , 
\\g £ b(-D)(A N , £ - Cl )- 1 - Fb(D)(A° N - Ciy^MO^MO) < C 7 p b (C)e 1/2 . 
The constants C7 and C7 depend only on the data m and ||Allx^oo • 


6.4. Estimates in a strictly interior subdomain. In a strictly interior 
subdomain O' of the domain O it is possible to obtain error estimates in 
H 1 of sharp order with respect to e. The following result was proved in 
Theorems 12.4 and 14.10 of |Su 8 l . 
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Theorem 6.5 ( }Su8j ). Suppose that the assumptions of Theorem 6.2 are 
satisfied. Let O' be a strictly interior subdomain of the domain O, and let 
5 = dist {O'; dOj. 

1°. Let ( € C \ M_|_ and |£| ^ 1. Then for 0 < e ^ £i we have 
|| {An,e C-0 1 ~ (S\° N - (I) 1 — eK N (e] 0\\ L2 (O)^ H i(O') 

^ (Csd- 1 +C 9 )c(<j)) 6 £, 

ii/6(d)(^, £ - (iy 1 -rs e b(p)Po(A% - ar'iiz^HMO') 

< (C^- 1 +C 9 )c{fife. 

The constants C 8j Cg, C 8; and Cg depend only on the initial data EH- 
2°. Now, let C € C \ [q,,oo), ( / 0. Denote p\,(C) := c (^)Pb(C) + 
c(i?) 5,/2 /9t > (C) 3 ^ 4 . For 0 < e ^ £\ we have 

II C^ >e - Cl)- 1 - (A% - a)’ 1 - eK N (e- C)V\\ L2 (0)^Hfi0') 

^ (CwS-'pfiO +C n c(#) 1/2 P fi0 5/4 ) s, 

\\g £ b(B)(A N , £ - Cl)' 1 - gFS e b(D)Po(A% - a)" 1 ^ \l 2 (0)^l 2 (0') 

^ pio<5 _ Vb(C) +Cnc(i?) 1/2 ft(C) 5/4 ) e. 

The constants C\g, C\ \, C\g, and C\ \ depend only on the data (J6.2D . 

The following result was proved in Theorems 12.5 and 14.11 of |Su8] : it 
concerns the case where A £ L^. 

Theorem 6.6 ( [Su8j ) . Suppose that the assumptions of Theorem 6.5 are sat¬ 
isfied. Suppose that the matrix-valued function A(x) satisfies Condition 2.5. 
Let K^{e\C) be the operator (16.101) . 

1°. Let C € C \ R+, |C| ^ I- For 0 < e ^ £\ we have 

\\(An,e ~ Cl) 1 - (A% - Cl) 1 ~ £K n( £ 'iC)\\l 2 (0)^h 1 (0') 

^ (C 8 <5 -1 + Cg)c((f>) 6 E, 

\\<fb(B){A N , e -Cl)' 1 -g £ b(D)(A° N -Cir^Wo^LW) 

^ (C 8 (5 _1 + Cg)c(fi) 6 E. 

The constants C 8 and C 8 are the same as in Theorem 6.5. The constants Cg 
and Cg depend only on the initial data (UP and IIAHi^ . 

2°. Let C € C \ [q,, oo) ; C A 0. For 0 < e ^ £\ we have 

\\{An,£~CI) 1 -(A° n -CI) 1 ~ £K %( £ ;C)\\l 2 (0)->-h 1 (0') 
^Cio<5 _1 a,(C) +Cnc(^) 1/2 ft(0 5/4 ) e, 

\\g £ b(D)(AN,E - Cl)- 1 -g £ b(mA% - Cir^Wo^m 
< {c w 8~ l pfiC) + Cnc(i?) 1/2 p b (C) 5/4 ) e. 










44 


YU. M. MESHKOVA AND T. A. SUSLINA 


The constants Cio and C 10 are the same as in Theorem 6.5. The constants 
C\\ and C\\ depend only on the data O) and HAHi^. 

§7. Homogenization of the operator exponential e~ AN ’ A 

7.1. The properties of the operator exponential. We start with the 
following simple statement about estimates for the operators e~ AN ' et V and 
e -A N t-p j n var j ous operator norms. 

Lemma 7.1. Suppose that the assumptions of Theorem 6.1 are satisfied. 
Let V be the orthogonal projection of 7^(0; C n ) onto the subspace (15.81) . 
Then for t > 0 and e > 0 we have 


~ AN ' et V\\ L2 (p)^L 2 (p) ^ 

(7.1) 

- AN ’ et nL 2{ 0)^Hpo) ^ c - 1/2 t-V 2 e- c ^/ 2 , 

(7.2) 

~ A N t V\\ l 2 (0)^l 2 (0) < e~ cA , 

(7.3) 

- A ^n uo^HHO) < c-^t-^e-^' 2 , 

(7.4) 

ANt V\\ l 2 (0)^h 2 (0) ^ c° t 1 e Cbt//2 , 

(7.5) 


where c° = c°(l + 1 ). 

Proof. Estimates (17.11) and (17.31) follow directly from (15.111) . 

By (EJ3D, 

We-^VW l 2(0 )^i (0 ) < c; 1/2 \\A]iy AN ^V\\ L2{0) ^L 2 (oy (7.6) 

Similarly to (13.71) . 

\\A)ly AN ’ A V\\ L2{0 )^L 2{ 0) < sup M 1/2 e-^ < (7.7) 

Relations (17.61) and (17.71) imply (17.21) . Similarly, from the inequality 

||(.<) 1 / 2 e-^T|| Il(0) ^ M0) < r'/V * 1 ' 1 ( 7 . 8 ) 

and (15.111) we deduce (17.41) . 

Next, by (15.71) . 

IMjv) 1 'P\\l 2 (0 )^h 2 (0) ^ c ° sup(/r + l)/i 1 ^ c°(l + x ) = c°. 
Similarly to (13.91) . this yields 

\\ e ~ ANt 'P\\L 2 { 0 )^H 2 ( 0 ) ^ c°\\A%e~ ANt 'P\\ L 2 ( K 0 )^L 2 ( 0 ) ^C°t 1 e _Cb ^ 2 , 
which proves (17.51) . □ 









HOMOGENIZATION OF INITIAL BOUNDARY VALUE PROBLEMS 


45 


7.2. Approximation of the operator e An .P j n the L2(0;C n )~ 
operator norm. Now, we prove the following result. 

Theorem 7.2. Suppose that the assumptions of Theorem 6.1 are satisfied. 
Then for 0 < e ^ E\ we have 

||g-Av,et _ e- A °^\\ L2{0) ^ L2i0) ^ C l 2 e(t + £ 2 )~ 1 / 2 e~ Cbt/2 , t > 0. (7.9) 

The constant C \ 2 depends only on the data (16.211 . 

Proof. Since I — V is the orthogonal projection of L 2 (0\ C n ) onto Z, by 
(15.101) . we have 

(j -V) = e~ A n\l -V) = I-V. 

Hence, 

e -A N ,et _ e -A%t _ ^ e -A N , e t _ e -A%t' J v (7. 10 ) 

Together with (15.111) this implies 

e -A N ,et_ e -A%t = _A_J~ e -tt ((^ _ c/) -l _ {A 0 n _ dC (7.H) 

Here 7 C C is the positively oriented contour consisting of two rays: 

7 = {C € C : Im £ ^ 0 , Re ( = Im ( + q,/ 2 } 
U{C€C:ImC<0, ReC = —Ini ( + q,/2}. 

Denote q, := max{l; y/5q,/2}. For ( € 7 and |£| ^ q, we apply estimate 
(16.411 . For ( £ 7 and |£| > q we use (16.311 . Similarly to the proof of Theorem 

3.2, this implies (17.91) . □ 


7.3. Approximation of the operator e An P in the (. L 2 —>• /7 1 )-norm. 

We introduce the operator 

JC N (t-e) = Ro[A e }S e b(B)Poe- A ^V. (7.12) 

The continuity of the operator /C^(t;e) from L 2 (0;C n ) to 17 1 (0; C n ) for 
t > 0 is checked similarly to the continuity of the operator (13.181) . 

Using Theorem 6.2, we obtain the following result. 

Theorem 7.3. Suppose that the assumptions of Theorem 6.2 are satisfied. 
Let /Cjv(i;e) be the operator (17.12[) . Then for 0 < e ^ E\ and t > 0 we have 

\\ e AN ' et ~ e ANt - £ ^n{P£)\\l 2 {0)^h 1 {0) 17 I3j 

+ et- 1 )e-*V 2 , ' " 

||/6(D)e-^^ - g £ S £ b(B)P 0 e- A ^V\\ L 2 { 0 )^ L 2 ( o) 

< Ciafc 1 / 2 *- 3 / 4 + er 1 )e-* t / 2 . 


The constants C 13 and C 13 depend only on the data 
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Proof. Similarly to (17.111) . we have 
e ~A N , e t _ e -A° N t _ e )C N (t- £ ) 

= e ~ (t - ar 1 - (A - a ) -1 - zk n {e- qv) d(. (7,15) 

Here K^{e\ (f) is the operator (16.51) . Note that, by (I5.10p . 

(A N ,e ~ Cl r 1 - (A - a )’ 1 - eK N {e- QV ^ 

= ((^, £ - Cl)- 1 - (A - a)" 1 - eK N {e- 0) V. 1,1,0 

For C € 7 and |C| q, we apply estimate (16.81) . For ( £ 7 and |C| > q, 
we use (16.61) and (|7.16l) . Similarly to the proof of Theorem 3.3 this implies 

(EH. 

Let us discuss the proof of estimate (17.141) . Since I — V is the orthogonal 
projection of onto Z = Ker 6 (D) = Ker^4jv,e, then 

g £ b{D)e~ AN A = g e b{ D)e~ Arf •'*?>. (7.17) 


Hence, 

g £ b(D)e~ AN ’V - 5 £ 5 £ 6(D)P 0 e'' 4 ^ t P 

= {9 £ b(P)(Me - Cl)- 1 -g £ S e bCD)P o (A 0 N - CI)- l )VdQ 

(7.18) 


Besides, we have 

g £ b(D)(A N ,e ~ (IQ 1 - T/S £ b(B)Po(A% - CQ^V ^ 

= (g £ b(D)(A N , £ - ClA ~ g e S e b(B)Po(A 0 N - CL) -1 ) V. ' ''' 

For C € 7 and |C| q, we apply (16.91) and (|7.19l) . For C £ 7 and |C| > q, we 
rely on (16.7|> . This leads to (17.141) . □ 

Remark 7.4. Taking Remark 6.3 into account, it is possible to prove esti¬ 
mates of the form (17.91) . (17.131) . and (17.141) with e -Cl,i / 2 replaced by 
where k > 0 is arbitrarily small number. Indeed, for sufficiently small e one 
can transform the contour of integration so that it will intersect the real axis 
at the point g® ~ Herewith, all constants in estimates of the form (|7.9j) . 
EH, and (17.141) will depend on k. 

7.4. Estimates for small t. Note that for 0 < t < e 2 instead of Theorem 
7.3 it is better to use the following simple statement (which is valid, however, 
for all t > 0 ). 
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Proposition 7.5. Suppose that the assumptions of Theorem 6.1 are satis¬ 
fied. Then for t > 0 and e > 0 we have 


lie AN ^-e A ° Nt \\L 2 (0)^m(0)<Cut 1/2 e ° A/2 

, (7.20) 

\\g £ b(D)e~ AN ’ et \\l 2 (0 )->l 2 (0) ^ Cut^ 2 e~ Ci,t ^ 2 , 

(7.21) 

|| 5 0 6(D)e-^‘|| L2(o) ^ i2(o) ^ Cut-^e-^ 2 . 

(7.22) 

Here C u = 2£ 1 _1/2 ||£T 1 ||^ and C u = WqWl* ■ 


Proof. Estimate (|7.2(jp follows from (|7.2p. (|7.4H. (17.101). and 
the constant q,. 

expression for 

By dzm 



\\g £ b(D)e An ^\\l 2 (o)^l 2 (o) < llslli,^ An ^V\\l 2 {0)->l 2 (0) ■ 

Together with (17.711 this implies (17.2111 . In a similar way, we check estimate 
(I7.22p . taking (|1.14l) and (17.811 into account. □ 

7.5. The case where A € Loo. Suppose that the matrix-valued function 
A(x) satisfies Condition 2.5. Then the operator 

JC° N (t-,e) = [A £ ]b(D)e~ A ^ (7.23) 

is continuous from L 2 ((D-,C n ) to 77 1 (0;C n ) for t > 0. 

Theorem 7.6. Suppose that the assumptions of Theorem 6.1 are satisfied. 
Suppose that Condition 2.5 is satisfied. Let Kfi N {t]£) he the operator (17.23p . 
Let g be the matrix-valued function 0 . Then for 0 < e ^ E\ and t > 0 
we have 

\\e~ AN < et - - £)C 0 N (t ; s)\\ L2 (0)^m(0) < Ci 5 (e 1 / 2 t ” 3/4 + £t~ l )e~ cA/2 , 

||/ 6 (D)e-' 4 ^* - g £ b(D)e- A ^\\ L2{0) ^ L2 (0) ^ C 15 (e l / 2 t ~ 3 / 4 + gf 1 )^/ 2 . 

The constants C 15 and C 15 depend only on the data dS 2 D and HAHi^ . 

Proof. The proof is based on application of Theorem 6.4 and repeats the 
proof of Theorem 7.3 with some simplifications. □ 

7.6. The case of the smooth boundary. It is also possible to remove the 
smoothing operator S e in the corrector under the assumption of additional 
smoothness of the boundary. In this subsection we consider the case where 
d ^ 3, because for d ^ 2 Theorem 7.6 is applicable (see Proposition 2.7(1°)). 

Lemma 7.7. Suppose that k 2 is integer. Let O C R d be a bounded 
domain with the boundary 80 of class C* -1 ’ 1 . Then for t > 0 the operator 
e -A° N t-p con ti nuous from L2(0]C n ) to H s (0;C n ), 0 < s ^ k, and the 
following estimate is valid: 

j| e - A N t 'p\\ l 2 (o)^h s (o) ^ S//2 e Cbt//2 , t > 0 . 


(7.24) 
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The constant'la’s depends only on s, ao, on, Uglily, ||g 1 ||l 00 ; the constants 
Ci and Co from the inequality 4E2D, the constant Ci from (HD, and the 
domain O. 

Proof. The proof is similar to that of Lemma 3.7. It suffices to check estimate 
(17.241) for integer s € [0, k]. For s = 0,1,2 estimate (17.241) has been already 
proved (see Lemma 7.1). 

So, let s be integer and let 2 ^ s ^ k. By [McL, Chapter 4], we conclude 
that the operator {A° N )~ S / 2 V is continuous from 7,2(0; C n ) to H s (0; C n ). 
Herewith, 

\\{A° n ) s / 2 V\\l 2 (o)^hs(o) (7-25) 

where the constant ffg depends only on s, ao, oq, H^Hloo, \\d _1 ||l 00 , the 
constants Ci and Co from the inequality (15.21) . the constant Ci from 
(15.91) . and the domain 0. Similarly to (13.321) . this implies (|7.24p with 
= %{s/e) s / 2 . □ 

Using the properties of the matrix-valued function A(x) and the operator 
S £ , it is possible to deduce the estimate for the difference of the correc¬ 
tors (|7.12|) and (17.231) from Lemma 7.7 under the assumption of additional 
smoothness of the boundary. 

Lemma 7.8. Let d ^ 3. Suppose that 0 C is a bounded domain with the 
boundary 80 of class C d Z 2,1 if d is even and of class C l ( d + 1 )/ 2 ’ 1 if d is odd. 
Let /Civ(t;£) be the operator (17.121) . and let /C^(t;e) be the operator ()7.23D . 
Then 

|| ALjv (i; e) — JC%(t.; £)||l 2 (o)-)-h1 (C)) ^ ^d{t d ^ 1 ^ 2 )e C|>t//2 , f > 0. 

(7.26) 

The constant ’rfd depends only on d, m, ao, a±, HlU 1 IIloo? the con¬ 
stants Ci and C 2 from the inequality (15.21) . the estimate Ci from (15.91) . and 

the domain 0. 


Lemma 7.8 and Theorem 7.3 imply the following result. 


Theorem 7.9. Suppose that the assumptions of Theorem 6.1 are satisfied, 
and let d ^ 3. Suppose that the domain 0 satisfies the assumptions of 
Lemma 7.8. Let KP N (t\e) be the operator (|7.23l) . Let g be the matrix-valued 
function dia. Then for 0 < e ^ £1 and t > 0 we have 


e A N>e t _ e A N t _ e]C 0 N (t-,£)\\ L2 ( O ^ H i( O ) 

^ ^*(£ 1/2 f -3/4 + £t _1 + st~ d/4 ~ 1/2 )e~ C],t/2 , 


Ub(D)e~ A ^* - 9 £ h{D)e- A >\\ L2{0) ^ L2{0) 
^ ^(£ 1/2 t -3 / 4 + £t -1 + £t _rf / 4—4 / 2 )e —Cb */ 2 . 
The constants and depend only on the data (16.21) . 


(7.28) 
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The proofs of Lemma 7.8 and Theorem 7.9 are postponed until Appendix 
(see §9). As in the case of the Dirichlet condition, it is convenient to apply 
Theorem 7.9, if t is separated from zero. 

Remark 7.10. In Lemma 7.8, instead of the smoothness condition on dO, 
one could impose an implicit condition on the domain: suppose that a 
bounded domain O with Lipschitz boundary is such that estimate ()7.25j) 
with s = d/2 + 1 is valid. For such domain all the assertions of Lemma 7.8 
and Theorem 7.9 are true. 

7.7. Special cases. Special cases are distinguished by the arguments sim¬ 
ilar to the arguments of Subsection 3.7. Theorems 7.3, 7.6 and Propositions 
1.2, 1.3, 2.7(3°) imply the following result. 

Proposition 7.11. Suppose that the assumptions of Theorem 7.3 are sat¬ 
isfied. 

1°. Let g° = ~g, i. e., relations ([1.121) are satisfied. Then /C 7 v(t;e) = 0, and 
for 0 < e ^ Ei we have 

\\e~ AN ’ et ~ e^||L a (0)-nP(0) < Ci 3 (e 1/2 f 3/4 + ef^e"^/ 2 , t > 0. 

2°. Let g° = g, i. e., relations (11.131) are satisfied. Then for 0 < e ^ £i and 
t > 0 we have 

||/6(D) e -^ - g 0 b(D)e- A ^\\ L2(O) ^ L2 (o) < C 15 (£ 1/2 f -3/4 + et- l )e- c ^ 2 . 

7.8. Estimates in a strictly interior subdomain. In a strictly interior 
subdomain O' of the domain O , using Theorem 6.5, we obtain estimates of 
sharp order 0 {e) (for a fixed t > 0 ). 

Theorem 7.12. Suppose that the assumptions of Theorem 7.3 are satis¬ 
fied. Let O' be a strictly interior subdomain of the domain O, and let 
5 = dist {O'; 80}. Then for 0 < £ ^ E\ and t > 0 we have 

|| e A N , s t _ e A N t _ eK.N(t;£)\\ L 2 ^ 0 ^ H i^ 0 ^ (729) 

^ (Ci6<5 -1 +Ci7)ef 1 e -Cb */ 2 , 

\\fb{B)e- AN ^ - FSeb^Poe-^Vh^L^ 

^ (C 16 r 1 +Ci 7 )et- 1 e- c ^/ 2 . 

The constants C\q, Cyj, C\&, and Cn depend only on the data (16.21) . 

Proof. The proof is similar to that of Theorem 7.3 and relies on Theorem 
6.5 and relations (17.151) . (17.161) . (17.181) . and ([7.191) . □ 


In the case where A E L^, Theorem 6.6 implies the following result. 

Theorem 7.13. Suppose that the assumptions of Theorem 7.12 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Let 
















50 


YU. M. MESHKOVA AND T. A. SUSLINA 


/C^-(i;e) be the operator (17.231) . Then for 0 < e ^ E\ and t > 0 we have 
|| e A N ,et _ e ■ A N t — e)C%(t-,s)\\ L 2 (o)^Hi( 0 ')^(C i6<5 1 +Cn)et l e Cb ^ 2 , 

\\fb(D)e- AN ^ - g e biP)e- A >\\ L2{0) ^ L2{0l) < (C^cT 1 + C^eT 1 e ^/ 2 . 

The constants C±q and C\o are the same as in Theorem 7.12. The constants 
C i7 and C\i depend only on the data m and 

It is possible to remove the smoothing operator S £ from the corrector 
in estimates of Theorem 7.12 without Condition 2.5. We consider the case 
d ^ 3 (in the opposite case, by Proposition 2.7(1°), Theorem 7.13 is ap¬ 
plicable). For t > 0 the operator e~ A ° Nt V is continuous from 7^(0; C n ) to 
77 2 (0; C n ), and estimate (I7.5P holds. Moreover, the following property of 
,,interior regularity” is true: for t > 0 the operator e~ A n tr P is continuous 
from 7^(0; C n ) to 77 s (0'; C n ) for any integer s ^ 3. We have 

||e _ • A ° Nt V\\ L 2 { 0 )^Hs( 0 >) < <t _1/2 (<T 2 + t > 0. (7.30) 

The constant depends on s, d, ao, cni, H^Hloo! I| 9 _1 ||l 0 o> the constants 
Ci and £2 from the inequality (15.21) . the constant £1 from (15.91) . and the 
domain 0. 

The way of the proof of estimates (17.30[) is the same as in the case of the 
operator A ° D ; see comments in Subsection 3.8. The difference is that instead 
of Lemma 3.1 one should use Lemma 7.1. 

The following statement is deduced from (17.301) by using the properties 
of the matrix-valued function A(x) and the operator S £ . 

Lemma 7.14. Suppose that the assumptions of Theorem 7.12 are satisfied, 
and let d ^ 3. Let /C° r (t; e) be the operator (17.231) . Let r\ be defined by (11.11) . 
Then for 0 < e ^ (4?’i) _1 5 and t > 0 we have 

||£jv( 7;e) — £jv(7; e )\\ L 2(p)^Hi(p') ^ ^d(t 1 +t ^ 2 (<5 2 + t l ) d ^ A )e c A/ 2 m 

(7.31) 

The constant’io'f depends only on d, ao, a\, 115 1 ||l 00 ; the constants 

£1 and £2 from the inequality m, the constant £1 from m and the 
domain 0. 

Lemma 7.14 and Theorem 7.12 imply the following result. 

Theorem 7.15. Suppose that the assumptions of Theorem 7.12 are satisfied, 
and let d ^ 3. Let KP n {L,e) be the operator (17.231) . Then for 0 < e ^ 
rninjei; (4ri) _1 <5} and t > 0 we have 

|| e A N , e t _ e A N t __ er/C^- (t; £)\\l 2 (0)^h 1 (0') ^ ^d^hdi^', t)e Cbt//2 , (7.32) 

\\g e b(D)e~ AN ^ - rb(D)e- A ° Nt \\ L 2 ( 0 )^L 2 ( 0 ') < V d eh d (5- t)e ~ c ^ 2 . (7.33) 

Here hd( 6 ;t) is defined by (13.41 [) . The constants % and % depend on the 
data (16.21) . 
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The proofs of Lemma 7.14 and Theorem 7.15 are postponed until Appen¬ 
dix (see §9). It is convenient to apply Theorem 7.15 if t is separated from 
zero. 

§8. Homogenization of the solutions of the second initial 

BOUNDARY VALUE PROBLEM 

In this section, we apply the results of §7 to homogenization of the so¬ 
lutions of the second initial boundary value problem for parabolic systems 
with rapidly oscillating coefficients. 

8.1. The problem for the homogeneous parabolic equation. We 

start with the case of homogeneous equation and nonhomogeneous initial 
condition. We are interested in the behaviour for small e of the generalized 
solution of the problem 

dUe 0 ?’ t:) = -6(D)V(x)6(D)u £ (x,f), x € O, t> 0, 

< u £ (x,0) = ip(x), xeO, (8.1) 

d^u £ (-,t)\do = 0, t> 0. 

Here (f € 7/2(0; C n ). The symbol df stands for the conormal derivative 
corresponding to the operator 6(D)*g £ (x)6(D). The formal differential ex¬ 
pression for df is given by d £ u(x) := 6(iz(x))*g £ (x)6(V)u(x). 

The solution of problem (18.11) is represented as u £ (-,f) = e^^ N ’ st ip. Thus, 
the question about homogenization of solutions is reduced to the question 
about approximation of the operator e~^ N P studied in §7. 

The corresponding effective problem has the form 

= -K D )V&(D)u 0 (x,t), xeO, t > 0, 

< u 0 (x,0) = <p(x), x<E0, (8.2) 

^u 0 (-,t)|ac> = 0, t> 0. 

We have u £ (-,t) = e _-/liv ’ e4 <£>, uo(-,t) = Therefore, Theorems 7.2 

and 7.3 together with the identity e^^^Vip = 'Pe~- Ab °N t Lp = Vvlq(-, t) imply 
the following result. 

Theorem 8.1. Suppose that the assumptions of Theorem 6.2 are satisfied. 
Let u £ and uo be the solutions of problems (El]) and (18.21) . respectively, 
where tp G 7/2(0; C"). Then for 0 < e ^ E\ and t^ 0 we have 

||u £ (-,t) - u 0 (-,f)|| L 2 ( e>) < C 12 e(t + e 2 ) _1/2 e“ Cbt/ -||(p|| L2(0) . 

Denote uo := PoVuq, p £ := g £ b( D)u £ . Then for 0 < £ ^ E\ and t > 0 we 
have 

||u £ (-,f) - u 0 (-,f) - eA £ S £ 6(D)u 0 (-, ^11^1(0) 

IIPeM) -g E S £ b(B)u 0 (-,t)\\ L ^ o] ^ C 13 (e 1/2 t~ 3/4 + et _1 )e _Cbt/2 ||^|| L2 ( 0) . 






52 


YU. M. MESHKOVA AND T. A. SUSLINA 


In the case where A £ Theorem 7.6 leads to the following result. 

Theorem 8.2. Suppose that the assumptions of Theorem 8.1 are satisfied. 
Suppose that Condition 2.5 is satisfied. Then for 0 < e ^ E\ and t > 0 we 
have 

||u e (-,f) -uo(-,t) -eA £ &(D)u 0 (-,t)||tfi (c >) 

< Ci5(e 1/2 i _3/4 + et _1 )e _Cbt/2 ||^|| L2(0) , 

IIPeM) -5 £ &(D)uo(-T)|| l 2 (0) ^ C 15 (e 1/2 t~ 3/4 + et _1 )e _c ^ /2 ||^|| l 2 (0)- 

In the case of sufficiently smooth boundary from Theorem 7.9 we obtain 
the following result. 

Theorem 8.3. Suppose that the assumptions of Theorem 8.1 are satisfied, 
and let dfi 3. Suppose that 80 € C d Z 2,1 if d is even and dO € C^ d+1 ^ 2,1 if 
d is odd. Then for 0 < £ ^ and t > 0 we have 

||ue(-,t) -u 0 (-,t) -eA £ 6(D)uo(-,t)|| H i (0) 

< ^(e 1/2 t~ 3/4 + ef _1 +ei _d/4 " 1/2 )e _Cbt/2 ||v?|| L2 ( 0 ), 

\\p £ (-,t) -g £ b{B)u 0 (-,t)\\ L2(O ) 

< ^(e 1 / 2 f 3/4 + ef _1 + £t~ d/A ~ ll2 )e~ c ^ tl2 \\if\\ L2{0) . 

Special cases are distinguished with the help of Proposition 7.11. 

Proposition 8.4. Suppose that the assumptions of Theorem 8.1 are satis¬ 
fied. 

1°. If g° = g, i. e., relations (|1.12D are satisfied, then for 0 < s ^ E\ and 
t > 0 we have 

||u e (-,f) - uo(-,t)|| H i (0) ^C n (e 1,2 r 3/A + ^ _1 )e _Cbt/2 11 ^ 11 ^( 0 )■ 

2°. If g° = g, i. e., relations (11.131) are satisfied, then for 0 < e ^ E\ and 
t > 0 we have 

IIPeM) -Po(-,*)||l 2 (0) ^Ci5(e 1/2 f 3/4 + et _1 )e" c ^ /2 ||(p|| i2(0) , 
where po = g ,0 6(D)uo. 

Approximation of the solutions in a strictly interior subdomain is deduced 
from Theorem 7.12. 

Theorem 8.5. Suppose that the assumptions of Theorem 8.1 are satis¬ 
fied. Let O' be a strictly interior subdomain of the domain O, and let 
5 = dist {O'; dO}. Then for 0 < £ ^ £\ and t > 0 we have 

||u e (-,t) - u 0 (-,t) - eA £ S e 6(D)u 0 (-, 

^ (Ci6<5 l +Cn)et 1 e C|,t//2 ||<p|| l 2 (0)> 

l|Pe(-,*) -? £ 5’^(D)G 0 (-,i)||L 2 (o') < (Ci 6 5 _1 +Ci 7 )et _1 e _c ^ /2 ||^||L 2 (0)- 

In the case where A £ L a 0 Theorem 7.13 implies the following result. 
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Theorem 8 . 6 . Suppose that the assumptions of Theorem 8.5 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Then 
for 0 < e ^ £\ and t > 0 we have 

||u £ (-,t) - u 0 (-,f) - eA £ 6 (D)u 0 (-,t)||^i( O /) 

^ (Ci6<5 1 +Cn)et 1 e c '° t ^ 2 \\f\\L 2 {0)i 

IIPeM) - <f&( D ) u o(u 011 ^ 2 ( 0 ') ^ (Ci6<J _1 +Cn)£t~ 1 e~ c '’ t/2 \\ip\\L 2 (py 

Finally, for dfi 3 Theorem 7.15 implies the following result. 

Theorem 8.7. Let d 3. Under the assumptions of Theorem 8.5 for 
0 < e ^ rnin{ei; (4rr) _1 <5} and t > 0 we have 

||ue(-,t) - uo(-,t) - eA £ 6(D)u 0 (-,t) 11^1(0/) < t &d£h d (5\ i)e -Ct,t/2 \\v\\ L2 ( 0 ) , 

||p e (-,t) - 5 £ 6 (D)uo(-,f)|| L2(c »/) < %eh (i (5;t)e _Cbt/2 ||^|| L2(cl) . 

77ere hfid]t) is defined by (I3.41|) . 

8.2. The problem for nonhomogeneous parabolic equation. Ap¬ 
proximation of the solutions in L 2 (0;C n ) and Sj p (T). Now, we con¬ 
sider the problem 

rXMx.f) = — 5(D)*g £ ( X )6(D)u £ (x, t) + F(x,f), x G 0, 0 < t < T, 

< u £ (x, 0) = f(x), x G 0, 
dfu £ (-,t)\oo = 0, 0 < t <T. 

(8.3) 

Here € L 2 (0;C n ) and F G Sj p (T) = L p ((0, T); L 2 (0; C n )), 0 < T < oo, 
for some 1 ^ p ^ oo. 

The corresponding effective problem is given by 

' dug(xj) _ _5(D)*gr 0 6(D) Uo (x, t) + F(x, t), x G 0, 0 < t < T, 

< u 0 (x, 0) = v?(x), x G 0, 

5°uo(-,i)bo = 0, 0 < t < T. 

(8 ' 4) 

We have u £ (-,t) = e~^ N > et f(-') + t) dt. Similarly, 

Uo(-,t) = e~- A N t f(-) + f*e~ A n( t ~ t )F(-,t)dt. By analogy with the proof 
of Theorem 4.8, we deduce the following result from Theorem 7.2. 

Theorem 8 . 8 . Suppose that the assumptions of Theorem 6.1 are satisfied. 
Letu £ be the solution of problem (|8.3D . where f G L 2 (0;C n ) and F G Sj p (T), 
0 < T ^ oo, with some 1 < p ^ oo. Let uo be the solution of the effective 
problem Then for 0 < e ^ £i and 0 < t < T we have 

l|u £ (-,i) —u 0 (-,i)|U 2 (0) ^Ci 2 e(t+£ 2 )- 1/2 e- c ^ /2 || ¥ 3|| i2{o) + c p 0(e,p)||F|| fip(t) . 

Here 9(e,p) is defined by (14.101) . The constant c p depends only on the data 
(16.21) and p. 
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The following result can be proved with the help of Theorem 7.2 by anal¬ 
ogy with the proof of Theorem 4.9. 

Theorem 8.9. Suppose that the assumptions of Theorem 6.1 are satisfied. 
Let u e be the solution of problem (18.31) . and let uo be the solution of problem 
(|8.4|) with ip € L2(0] C n ) and F € Sj p (T), 0 < T ^ oo, for some 1 ^ p < oo. 
Then for 0 < e ^ E\ we have 

K - u o||ij p (T) ^ V0(e,p)||^|| i2 ( O ) +Ci 8 e||F|| %(T ). 

Here @(e,p) is defined by (14.171) . and Cig = C\2(2tt) 1 I 2 c^ 1 ^ 2 . 

8.3. Approximation of the solutions in H l {0\ C n ). Now we obtain 
approximation of the solution u e of problem (18.31) in i7 1 (0; C n ). By analogy 
with the proof of Theorem 4.11, we deduce the following result from Theorem 
7.3 and Proposition 7.5. 

Theorem 8.10. Suppose that the assumptions of Theorem 6.2 are satisfied. 
Let u £ be the solution of problem (18.3jl . where ip £ L2(0]C n ) and F £ Sj p (T), 
0 < T ^ oo, for some 2 < p ^ oo. Let uo be the solution of the effective 
problem (18.41) . Assuming that t ^ e 2 , denote 

w e (-,t) :=e-< £2 u 0 (-,f-e 2 ) =e-<V(')+ £ e -A° N (t-t) F fj)dt. 

Jo 

(8.5) 

We put w e (-,t) := P 0 Vw £ (-,t), p £ = g £ b(D)u e . Then for 0 < e ^ E\ and 
e 2 ^t<T we have 

||u e (-,f) -uo(-,i) - eA £ S £ b(D)w £ (-,t)\\ H i {0) 

< 2 Ci3e 1/2 f 3/4 e _Cbt/2 ||<p||L 2 (o) + c P p(e,p)||F|| %(t) , 

\\Pe(-,t) -g £ S £ b(D)w £ (-,t )|| La ( 0) 

< 2 Ci3e 1/2 t _3/4 e _Cbi/2 ||<p||L 2 (o) + Cpp(e,p)||F|| %(t) . 

Here p(e,p) is defined by (14.241) . The constants l p and c p depend only on 
the data (62D and p. 

In the case where A £ L ^ Theorem 7.6 and Proposition 7.5 imply the 
following result similar to Theorem 4.12. 

Theorem 8.11. Suppose that the assumptions of Theorem 8.10 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. Then 
for 0 < £ ^ E\ and e 2 ^ t < T we have 

||u £ (-,i) -u 0 (-,t) -eA £ 6(D)w £ (-,t) 

< 2 <Ji5£ 1/2 t~ 3/4 e-^ t/2 \\ip\\ L2 (o) + c pP(^P)ll F IUp (t) , 

IIPeM) -<f&(D)w e (-,t)|| l 2(0) 

< 2 Ci 5 e 1 / 2 t _ 3 / V Cki/ 2 ||v?|| L2 ( 0 ) + c p p{£,p)\\F\\^ At) . 

The constants c p and c" depend only on the data (EH). ||A|| ioo , and p. 
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In the case of additional smoothness of the boundary Theorem 7.9 implies 
the following result (cf. Proposition 4.13). 

Proposition 8.12. Suppose that the assumptions of Theorem 8.10 are sat¬ 
isfied. Let d = 3 and p = oo. Suppose that dO € C 2,1 . Then for 0 < e ^ E\ 
and e 2 ^t < T we have 

||u £ (-,t) - uo(-,t) - eA e &(D)w £ (-,f)|| H i (0) 

< 2%V/2 t -3/4 + rf-5/ 4 )e- c ^ 2 M L2(0) + c , e 1 / 2 ||F|| floo(t) , 

\\Pe(-,t) -g e b(D)w e (-,t)\\ L2{0 ) 

< 2^|(e 1/2 f3/4 + £t -^)e - c ^ 2 || ¥ »|| La(0) + ?e 1 / 2 || F || flao(t) . 

The constants c 1 and c' depend only on the data 

Now we distinguish special cases. The following statement is proved with 
the help of Propositions 7.5, 7.11 and Theorem 8.10, by analogy with the 
proof of Proposition 4.14. 

Proposition 8.13. Suppose that the assumptions of Theorem 8.10 are sat¬ 
isfied. 

1°. If g° = g, i. e., relations () 1.12 f) are satisfied, then for 0 < e ^ E\ and 
s 2 ^t <T we have 

||u £ (-,f) - Uo(-,t)|| H l(o) 2Ci 3 e 1 /2 t -3/4 e - Cb t/2|| ¥ ,||^ (c)) 

+ c p p(£,p)||F(-,t)|| flpW . 

2°. If g° = g, i. e., relations (11.131) are satisfied, then for 0 < e ^ E\ and 
£ 2 ^t <T we have 

IIPeM) -PoM)IIl 2 ( 0 ) < 2Ci 5 e 1/2 t“ 3/4 e“ c ^ /2 ||^|| L2(0) + CpV(e,P)l|F|| flp(t ), 

where po = g°6(D)uo. The constant c p " depends only on the data (16.21) and 

p. 

8.4. Approximation of the solutions in a strictly interior subdo¬ 
main. Approximation of the solutions in a strictly interior subdomain can 
be derived from Theorem 7.12 (cf. Theorem 4.15). 

Theorem 8.14. Suppose that the assumptions of Theorem 8.10 are sat¬ 
isfied. Let O' be a strictly interior subdomain of the domain O, and let 
5 = dist {O'; dO}. Then for 0 < e ^ E\ and e 2 ^ t < T we have 

||u £ (-,f) -u 0 (-,t) -£A £ 5 e 6(D)w e (-, t)\\ H i {ol) 

< (Ci6<5 _1 +Ci 7 )et _1 e _Cbt/2 11^11^(0) + (k p( 5 _1 + kpcr(e,p)||F|| j = lp(t) , 

IIPeM) - g £ Seb(D)w £ (-,t)\\ L2 (0') 

< (CiecT 1 +C l7 )£t- l e- c '> t/2 \\ip\\ L2 ( 0) + (kpfT 1 + k^) C j(e,p)||F|| j j p(t) . 

Here cr(£,p) is defined by (14.39D . The constants k p , k p , k p , and k p depend 
only on the data (16.21) and p. 
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In the case where A G Loo from Theorem 7.13 we deduce the following 
result similar to Theorem 4.16. 

Theorem 8.15. Suppose that the assumptions of Theorem 8.14 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. For 
0 < e ^ S\ and e 2 ^ i < T we have 

||u £ (-,t) -uo(-,i) -eA £ 6(D)w £ (-,i)|| H i (0 /) 

< (Ci6<5 _1 +C 17 )er l e- C[ ° tl2 \\v\\ L2{0) + (kpcT 1 + kpcr(e,p)||F|| %(t) , 

\\Ps(-,t) -g £ b(D)w £ {-,t)\\ L2{0/) 

< (Ci6<5 _1 +Ci 7 )et- 1 e- c * t/2 \\(p\\ La ( 0) + (kpfT 1 + kpu(e,p)||F|| %(t) . 

The constants k p and k p are the same as in Theorem 8.14. The constants 
k p and \fi p depend only on the data ra, \\a\\ L oo , and p. 

8.5. Approximation of the solutions in ©p(T) = L p ((0, T); i7 1 (0; C n )). 
By analogy with the proof of Theorem 4.18, it is easy to check the following 
result. 

Theorem 8.16. Suppose that the assumptions of Theorem 6.2 are satisfied. 
Let u £ and uo be the solutions of problems (|8.3I) and (|8.4I) , respectively, 
where f G L2(0;C n ) and F € S) p {T), 0 < T ^ oo. Suppose that for 
e 2 ^ t < T the function w e (-,t) is defined by (18.51) . For 0 ^ t < e 2 we put 
w £ (-,t) = 0. Denote w £ := P 0 Vw £ (-,t). Let p £ = g £ b( D)u e . 

1°. Let 1 ^ p < 2. Then for 0 < e ^ E\ we have 

l|u £ - Uo - eA £ S' £ 5(D)w £ || l5p(r) ^ K,' p a(£,p)\\f\\ L2{0] + Ci 9 e 1/2 ||F|| %(T) , 

11 Per — 9 e S £ b(F>)y^e\\ Lp ((p, T )-, L 2 { 0 )) ^ K p a ( £ ) p) II T\\l 2 {0) +Ci 9 e 1/,2 ||F|| sj p (T)- 

Here a(e,p) is defined by (14.411) . The constants C\g and C\g depend only 
on the data (EH). The constants K p and k" depend on the same parameters 
and on p. 

2°. Let f = 0, and let 1 ^ p ^ oo. Then for 0 < e ^ £\ we have 

||u £ - u 0 - eA £ 5 £ 6(D)w £ || 0p(T) sC Ci 9 e 1/2 ||F|| %(T) , 

IIPe -? e 5' £ 6(D)w £ || I/p (( 0iT ). L2 ( O )) ^ Ci 9 e 1 ' /2 ||F|| i5p ( T ). 

Under the additional assumption that A G L^ we obtain the following 
theorem similar to Theorem 4.19. 

Theorem 8.17. Suppose that the assumptions of Theorem 8.16 and Con¬ 
dition 2.5 are satisfied. 

1°. Let 1 ^ p < 2. Then for 0 < e ^ s± we have 

l|u £ - Uo - eA £ 6(D)w £ ||g p(T) < t p a(£,p)\\f\\ L2{0) + C 20 e 1/2 ||F||f lp(T) , 

IIPe — S ,e K D ) w £llL p ((0,r);L 2 (O)) ^ £p a ( £ iP)\\ l f\\L2{0) +C 2 oe 1//2 ||F|| j = lp ( T ). 
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The constants C 20 and C 20 depend only on the data (|6.2I) ||A|| Loo . The 

constants t p and t p depend on the same parameters and on p. 

2°. Let = 0, and let 1 ^ p ^ 00 . Then for 0 < e ^ £\ we have 

||u e - 110 - eA £ 6 (D)w £ || (5p(T) ^ C 20 £ 1 / 2 ||F|| i -, p(T) , 

l|Pe - <f &(D)w £ || ip (( 0 ) x) ; L 2 (o)) < C 2 oe 1 / 2 ||F|| flp ( T ). 

The following result is similar to Proposition 4.20 and is deduced from 
Theorem 7.9. 

Proposition 8.18. Suppose that the assumptions of Theorem 8.16 are sat¬ 
isfied, and let d = 3. Suppose that dO E C 2,1 . 

1°. Let 1 ^ p < 4/3. Then for 0 < £ ^ £1 we have 

IIUf - Uo - eA £ 6 (D)w £ || ep(r) < tp£ 2A ^ 3/2 ||^|| L2(0) + C 2 i£ 1 / 2 ||F|| %(T) , 

llPe - 5‘' 6 (F))w £ || ip ((o,T);L 2 (e>)) ^ t p e 2 ^ p ^ ||<P||l 2 ( 0 ) +C 2 i£ 1 '-||F|| %(T) . 

The constants C 21 and C 21 depend only on the data m- The constants t' p 

and depend on the same parameters and p. 

2°. Let cp = 0, and let 1 ^ p ^ 00 . Then for 0 < £ ^ £-{ we have 

Il u e — u o - £A e 6 (D)w £ ||< 3 p ( T ) ^ Caie^HFlU^r), 

llPe - 7 eb (F>)w £ || Lp (( 0iT ) ;L2 ( O) ) < C 2 ie 1 // 2 ||F|| flp ( T ). 

8.6. Approximation of the solutions in L p ((0, T); H l (0'\ C n )). The 
following result is similar to Theorem 4.21 and can be checked by using 
Theorem 7.12 and Proposition 7.5. 

Theorem 8.19. Suppose that the assumptions of Theorem 8.16 are sat¬ 
isfied. Let O' he a strictly interior subdomain of the domain O, and let 
5 = dist{ 0 '; 00 }. 

1°. Let 1 ^ p < 2. Then for 0 < £ ^ £1 we have 
IIUe — U 0 - £A e S , £ 6 (D)w £ || ip (( 0 iT ).^l( C ,/)) 

^ (k p>5 1 + kp,)r(£,p)||<p|| i2 (£)) + (C 22 <5 1 + C 22 )£(| In e| + l)||F|| fip ( T ), 
llPe _ 5 o 5 e b(D)w e || ip (( 0iT ) ;i2 ( C )/)) 

^ (kp/J 1 + kp/)r(£,p)||y>||i, 2 (o) + {T 22 5 1 + ^ 22 ) e (l l ne l + l)l|F||^ p (T). 
Here r(e,p ) is defined by (14.541) . 

2°. Let ip = 0, and let 1 ^ p ^ 00 . Then for 0 < £ ^ £1 we have 
||u £ — u 0 — £A"S £ 6 (D)w £ || ip (( 0 !T ) ;i p( C >/)) 

^ (C 22 ^ 1 +C 22 )£(| lne| + l)||F|U p(T) , 
llPe - ff e 5' £ 6(D)w £ || Lp (( 0iT ) ;i2 ( O /)) < (C 22 <5 1 + C 22 MI lne| + l)||F|| fip ( T ). 

The constants C 22 , C 22 , C 22 , and C 22 depend only on the data (I6.2p . The 
constants k p /, k pl , k p /, and k!, are the same as in Theorem 8.14 (with p 
replaced hyp'). 
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Finally, in the case where A € Loo Theorem 7.13 and Proposition 7.5 
imply the following result similar to Theorem 4.22. 

Theorem 8.20. Suppose that the assumptions of Theorem 8.19 are satisfied. 
Suppose that the matrix-valued function A(x) satisfies Condition 2.5. 

1°. Let 1 ^ p < 2. Then for 0 < s ^ £\ we have 

Il u e - u 0 — eA^D^sll^^o^jL^o')) 

^ (k p /5 1 + kp/)r(e,p) 11^11^(0) + 1 + C 22 MI l ne l + l)l|F||j} p (TL 

IIPe ~ 9 £ b(D)'W £ \\L p ((0,T);L 2 (O / )) 

^ (kp/5 1 + kp,)r(e,p)||< pIIl 2 (0) + (C 22 5 1 + C 22 MI l 116 ! + 1)l|F||.$ p (x)■ 
2°. Let ip = 0, and let 1 ^ p ^ 00 . Then for 0 < e ^ £1 we have 

Il u e — u o — eA e &(D)w e || ip (( 0 iT ).#i( O /)) 

^ (<^22^ 1 + 1 )II F II%(T)! 

l|Pe ~ fif"&(D)w e || ip (( 0) T) ; L2(O')) ^ (^22^ 1 +^2) e (|l ne l + 1) ll F ll%(T) • 

The constants C 22 and C 22 are the same as in Theorem 8.19. The constants 
C 22 and C 22 depend only on the initial data (16.21) and HAjlx^. The constants 
k p i, kp,, k p i, and \s! pt are the same as in Theorem 8.15 (with p replaced by 

p'). 


Appendix 

§9. Removal of the smoothing operator 

In §9, we assume that d 3 and prove the statements about removing 
the smoothing operator S e in the corrector in the case of sufficiently smooth 
boundary (Lemma 3.8 and Theorem 3.9 in the case of the Dirichlet condition, 
Lemma 7.8 and Theorem 7.9 in the case of the Neumann condition) and in 
a strictly interior subdomain (Lemma 3.15 and Theorem 3.16 in the case of 
the Dirichlet condition, Lemma 7.14 and Theorem 7.15 in the case of the 
Neumann condition). 

9.1. The multiplicative property of the matrix-valued function 

A(x). 

Lemma 9.1. Suppose that the matrix-valued function A(x) is the T-periodic 
solution of problem Let dfi 3, and let l = d/2. 

1°. For 0 < e ^ 1 and for any v € T7 i_1 (M rf ; C m ) we have A £ v € L 2 (M rf ; C n ), 
and 

ll^° V llL 2 (R d ) ^ C® ll v lltff-i(Rd)- (9-1) 

2°. For 0 < e ^ 1 and for any v € C m ) we have A £ v g i7 1 (M d ; C n ), 

and 


ll A£ v||tfl( K d) < C [1) £ X ||V 


(9.2) 
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The constants C^\ C^ l \ and C^ depend on m, d, ao, ot\, H^Hloo? 
|| 5 - 1 || Loo; and the parameters of the lattice T. 


Proof. It suffices to check (19.111 and (19.21) for v € Co°(M d ; C m ). Substituting 
x = ey, e d// 2 v(x) = V(y), we have 

H A£v llL 2 (R d ) ^ [ l A ( e_lx )l 2 l v ( x )l 2dx = f l A (y)l 2 l v (y)l 2 ^y 

JR d J R d 

= E / |A(y)l 2 |V(y)l 2 rfy<E ll A l| 2 2 r(n)ll V l| 2 27 .,(n+a)J 

a gr •'° +a aer 

where r _1 + (r 7 )^ 1 = 1. The number r is chosen so that the continuous 
embedding C L 2r (Q) is valid, i. e., r = d(d — 2) _1 . Then 

ll A lll 2r (n) ^ c o||^|||fi(n); (9-4) 

where the constant cq depends only on the dimension d and the lattice T. 
For r chosen above we have 2 r' = d. From the continuity of the embedding 
H 1-1 (Q) C Ldifd) it follows that 

l|V|lL ( n + a) ^ c nll v ll^-i(n+ a )> ( 9 - 5 ) 

where the constant Cq depends only on the dimension d and the lattice T. 
Now, relations (|9.3D - (l9.5p imply that 

J Rd l A ( e_lx )l 2 l v ( x )l 2dx < c s^ c si IIA|llri ( S7) IIVll^z—i (]R d } - (9.6) 

Obviously, for 0 < e ^ 1 we have ||V11i( ir^) ^ Combining 

this with (19.611 and (11.181) . we obtain 


/ |A(e 1 x)| 2 |v(x)| 2 dx ^ cnc , n M 2 ||v||^ ! _ 1 

J R d 

which proves (19.11) with = cqCqM 2 . 
Next, by Lemma 1.5, 

l|D(A £ v)|| 2 2(Rti ) 


€ C 0 °°(ir;C m ), (9.7) 


^ 2e 


J^ d |(DA) £ (x)v(x)| 2 dx + 2 |A £ (x)| 2 |Dv(x)| 2 dx ^ 

2/?ie -2 [ |v(x)| 2 dx + 2(1 + j3 2 ) f |A £ (x)| 2 |Dv(x)| 2 dx. 

J R d JR d 

From (19.71) (with v replaced by the derivatives djv) it follows that 

^J A£ (x)| 2 |Dv(x)| 2 dx ^ cnc^Af 2 ||v||^ i(]Kd) , v£C 0 “(I rf ;n. (9.9) 

As a result, relations (19.7D (19.91) imply (19.21) with the constants C^ 1 ) = 
(2ft) 1 / 2 and C ^ = M{ 3 + 2 /3 2 ) 1/2 (cnc' n ) 1/2 . □ 


Using the extension operator Pq satisfying estimates (12.5p . we deduce the 
following statement from Lemma 9.1(1°). 
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Corollary 9.2. Suppose that the assumptions of Lemma 9.1 are satisfied. 
Then the operator [A £ ] of multiplication by the matrix-valued function A £ (x) 
is continuous from to L 2 ( 0 \C n ), and 

\\[^ £ ]\\h 1 - 1 {0)^L 2 (0) ^ (C (0) ) 1/2 cg 

9.2. Proof of Lemma 3.8. Suppose that the assumptions of Lemma 
3.8 are satisfied. Let uo(-,f) = e~' 4 » , Y>(.), where ip € L2(0;C n ). We 
put Uq(-, f) = Poii 0 (-,t). According to (13.181) and (13.271) . K.£)(t;£)(p = 
A £ 5 e 6(D)uo and Kfi D (t\e)<p = A £ 6(D)uo. We have to estimate the following 
function 

\\K D {fie)v - K° D (t-e)v\\m {0) = ||A £ ((S £ - I)b( D)5 0 ) (;t)\\m(o) . 

<\\A £ ((S e -I)b(-D)u 0 )(.,t)\\ mm . 


By Lemma 9.1(2°), we have 

\\JC D (t-, e)p> - JC° D {t; e)(p\\ H i (0) 

^ cWe- 1 !! ((S e - /)6(D)uo) (•,i)|| L2(Rd ) (9.11) 

+ C^\\ ((S £ - I)b( D)uo) l = d/2. 

The first term in the right-hand side of (19.111) is estimated by using 
Lemma 1.8 and (11.41) . (12.51) . (13.51) : 

e_1 || ((S £ - I)b{ D)uo) (-,t)||L 2 (Md) 

< ri||D6(D)u 0 (-,i)|| ia(R d ) ^ riai /2 ||u 0 (-,t)|| H 2 (Rd) (9.12) 

< riai /2 C'g 2) ||uo(-,t)|| i/ 2 (c)) ^ C (3) t _1 e“ c * t/2 ||^|| L2(CI) , 

where C^ 

In order to estimate the second term in the right-hand side of (19.111) . we 
apply the inequality ||S , e ||if!( R d)_ K f/i( R d) ^ 1 and (fLTD : 

II ((-Se — d)bfD) U 0 ) (-, ^)||iY'(K d ) ^ 2 IIK D ) u o(-^)llH i (R [i ) 

V2„~ , (9- 13 ) 

By (|2.5j) and Lemma 3.7, 

11 u 0 (• 5 i )ll,ff i + 1 (R d ) ^ cg + II u 0 (•, II (c>) . . 

^C'g +1) Q +1 f-( z+1 )/ 2 e- c * t / 2 ||^|| L2(0) . 

Relations (19.131) and (19.141) imply that 

|| ((5 e - /)6(D)uo) < cWt-V+We-^ 2 Ml 2{0 ), (9-15) 


where = 2a^ 2 C'g +1 ^C/ + i. 
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As a result, relations (19.111) . (19.121) . and (|9.15p imply that 
||/C D (t; e)p> - JC° D (t ; e)<p\\ H i (0) 

^ (c^c^h- 1 + e~ c ^ 2 M L2(0 ), 

which proves (13.331) with the constant = maxjC'WC'®; C^C^}. This 
completes the proof of Lemma 3.8. 


9.3. Proof of Theorem 3.9. Inequality (13.341) follows directly from (13.191) 
and (13.331) . Herewith, C* d = C 13 + C^. 

Let us check (13.351) . From (13.341) and (11.21) . (11.51) it follows that 

||/6(D) - £ A'6(D) e -^‘) || WOKI , ( „ 

v ' (9.16) 

< ll9lli„(<i«i) 1/2 C;(e 1/2 *- 3/4 +et- 1 + 

According to (11.21) . 

/6(D) (eA e 6(D)e-“ 4 n^ 

d (9.17) 

= g £ (b(D)Afb(T>)e- A D t + e ^2 g^^bjD.Dje-^ 1 . 

i,j =1 

The norm of the second summand in the right-hand side of (19.171) is esti¬ 
mated with the help of (1.5), Corollary 9.2 and Lemma 3.7: 


y, g £ b i A £ b j D i D 


;e A ° Dt 


*,l=i 


L 2 (0)^L 2 {0) 


^ e ll5l|Loo a l^(C ,( ' 0 ^) 1 ^ 2 C , O ^ll e Ar,t \\ L 2 (0)->H i + 1 (e>) 
< 2 e -c.t /2 


(9.18) 


where l = d /2 and C® = ||5 ||l O o q; i^(C ,( ' 0 ' ) ) 1 // 2 C , o ^Q+i- 

Now, relations ()9.16|) — (19. 181) and (1.9) imply (13.35[) with the constant 
C* d = lbl|Loc(rfai) 1/2 Cd + C^). This completes the proof of Theorem 3.9. 


9.4. Proof of Lemma 7.8. We prove Lemma 7.8 by analogy with the 
proof of Lemma 3.8. Let ip € L-2(0;C n ), and let ip = Vip. We put 
uo(-,f) = e~' 4 < v t <£(-) and uo(-,f) = Pc>uo(-,t). According to (17.121) and 
(17.231) . K.N{t\e)ip = A £ S' £ 6 (D)uo and KP N {t\e)ip = A £ 6(D)uo- We have 
taken into account that b(D)e^ A °N t = b(D)e^ A °N t V. 

Similarly to (19.101) . we obtain 


IIe)^ — £ )IliY 1 (O) / l|A £ ((S £ — I)6(D)u 0 ) (•, t)|| H i( R d). 
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Further considerations are similar to (|9.11j) - (l9.15jk instead of (|3.5j) we use 


m, instead of Lemma 3.7 we apply Lemma 7.7. This yields 

II £ )P - £jU*; e)v>lltfi(©) 


< (c^cPt~ l + e- c ^ 2 ||v.|| L2(0) 


where cffl = r\a^ 2 C$c° and = 2aq //2 cg +1 ^ + i. We arrive at (I7.26P 
with the constant % = maxjC'WC'®; Lemma 7.8 is proved. 

9.5. Proof of Theorem 7.9. Arguments are similar to the proof of Theo¬ 
rem 3.9. Inequality (17.271) follows directly from (17.131) and (17.261) . Herewith, 
= ^-13 + %• 

Inequality (I7.28P is deduced from (17.271) by analogy with (19.161) (19.181) . 
In addition, we use the identity 


d d 


e Y g e b i A £ b j D i D j e~^° Nt = e Y cfbitfbjDiDje-^V (9.19) 



and estimate this term with the help of Corollary 9.2 and Lemma 


7.7. We arrive at (|7.28l) with the constant = IlfflUoo (^ a i) 1//2 ^i + 
||< 7 || ioo a 1 d(C(°)) 1 /2cg- 1 ^ +1 . This completes the proof of Theorem 7.9. 

9.6. One property of the operator S e . We proceed to estimates in 
a strictly interior subdomain and start with one simple property of the 
operator S e . 

Let O' be a strictly interior subdomain of the domain O, and let 6 = 
dist {O'; dO}. Denote 

O" = {x € O : dist {x; 80} > 5/2}, O'" = {x € O : dist {x; 80} > 5/ 4}. 

Lemma 9.3. Let S e be the operator (11.191) . Let r\ be defined by (11.11) . Then 
for 0 < e ^ (dri)^ 1 ^ and for any function v € H s (R d -,C m ), s E Z + , we 


have 


ll-S'wlUqo") ^ Il v 


H s (0"')- 


Proof. According to (11.191) 


2 



D“v(x — sz) dz 


(9.20) 



Since 0 < er\ ^ <5/4, then for x E O" and z E LI we have x — ez E O'". 
Then, changing the order of integration in (19.201) . we obtain 
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9.7. The cut-off function x( x )* We fix a smooth cut-off function x( x ) 
such that 

X € C 0 °°(M d ), o ^ x( x ) < 1; x( x ) = 1, X e O'- 
suppx C 0"\ |D“x(x)| < Ss5~ s , \a\ = s, Vs € N. 

The constants <j s depend only on d, s, and the domain O. 


Lemma 9.4. Suppose that the function x( x ) satisfies conditions (19.211) . Let 
k € Z + . 

1°. For any function v € C m ) we have 

k 

llx v lltf fc (R d ) ^ (9.22) 

s=0 

where the constant C/. J depends only on d, k, and the domain O. 

2°. For any function v € F[ k+1 (M . d ; C m ) we have 

/k -\-1 

IIX V ll_H-fc+l/2(Rd) ^ Cfc+1/2 ( ^ ^ +1 

/ k 

x <A T, |I V II H°{0") 

V=o 

where the constant C'^ 1 ^ /9 depends only on d, k, and the domain O. 

Proof. Inequality (19.221) follows from the Leibniz formula for derivatives of 
the product x v and estimates for derivatives of x (see (19.211) 1. To check 
(|9.23|) . in addition one should use the obvious inequality 

H w ||^fc+i/2( R d) ^ ll w lliL fc + 1 (R d ) ll w lll? fc (R' i )! w £ H k+1 (U d -,C m ). 

□ 


v ll H“(0") 
1/2 


1/2 


(9.23) 


9.8. Proof of Lemma 3.15. Suppose that the assumptions of Lemma 
3.15 are satisfied. Let uo(-,t) = 4 <£>(•), where ip € L2(0-,C n ). By (11.151) 

and (13.81) . 

II Du o(-^)IIl 2 (C>) ^ c 0 1//2 \\(-A°d) 1//2 e ADt \\L 2 (0)^L 2 (0)\\ l P\\L2{0) (Q24) 

< C 0 1/2t ~ 1/2 e~ c * t/2 \\<p\\ l 2 (0)- 

From (13.51) it follows that 

11 u o (• 5 t)\\H 2 (o) < cJ _1 e _C * t/2 || < P||L 2( o ) - (9.25) 

Let uo = PoUo- According to (|3.181) and ()3.27l) . /Co(f;£)<£> = A £ S £ b(D)uo 
and JCpft-, e)(p = A e 6(D)uo- 

We have to estimate the function 


\\^D(t;e)(p - K^,(t;£)(p\\ H i^ ol) ^ ||A £ x ((^ - I)6(D)u 0 ) (•, t)|| H i (Rd) . 

(9.26) 
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By Lemma 9.1(2°), 

l|A £ X {{S e — I)6(D)u 0 ) (•, £)||#i ( Rd) 

< C [1) e- 1 ||x {{Se - /)6(D)uq) (•, t) || L2(Rd) (9.27) 

+ C^\\x {{S £ - /)6(D)uo) (■, t)\\ Hl[Rd) , l = d/2. 

The first term in the right-hand side of (19.271) is estimated with the help 
of Lemma 1.8 and (11.41) . (12.51) . (19.25p . Similarly to (|9.12l) . we have 


£ - 1 ||x((5 £ -/)6(D)5 0 )(-,t)|| i2 


< e X || {{S £ - I)b{ D)u 0 ) (•,t)|| L2(R d ) 




(9.28) 


We proceed to the second term in the right-hand side of (19.271) . Obviously, 
||x((S , e -/)6(D)uo)(-,t)|| H «(R- ) < ||x(5'e&(D)G 0 ) (■,i)llff' ( R«') 

+ ||x^(D)u 0 (-,t)|| i? i( R d). 


(9.29) 


To estimate the second term in the right-hand side of (19.291) . we apply 
Lemma 9.4 and (11.2j) . (11,51) . In the case of integer l = d/2 (i. e., d is even) 
we have 

i 

||x&(D)u 0 (-,t)||jyi ( r 4) ^ C z (6) (dai) 1/2 ^(5' ( ' _s) ||Duo(-,t)||^ S (o")- (9-30) 

s =o 

In the case of half-integer l = d/2 = k + 1/2 (i. e., d is odd) we have 

/k-\-l \ 1/2 

\\ X b{n)u 0 {;t)\\H‘m ^cf/da,) 1 / 2 (J2^ k+1 -^\\Bu 0 {.,t)\\Hs { o")) 

^s =0 ' 

/ k \l/2 

{k T, |l Du o('! i )llH^(0") j 

v=o ' 

(9.31) 


The norms of Duo(-, t) in L2{0) and H l {0) have been estimated in (19.241) 
and (19.251) . By (13.371) (with O' replaced by O"), 

l|Du 0 (-, t)\\ h s { 0 ") ^ C' +1 2 s t 1//2 (i5 2 + t 1 ) S//2 e c * t//2 1|IIz. 2 (0)’ s ^ 2. 

(9.32) 

Using (I9.24p . (I9.25p . and (19.301) (19.321) . we arrive at 

||x6(D)u 0 (-,t)|| ff i( R d) ^ c , (7)r 1/2 (<r 2 + t _1 ) d/4 e _c * t / 2 |MIl 2 ( 0 )- (9-33) 

The constant C^ depends on d, qo, a i, ||5'||i ;oo , ||<7 _1 ||l cx ,> and the domain 
O. 

To estimate the first term in the right-hand side of (19.291) . we apply 
Lemmas 9.4 and 9.3. Assume that 0 < £ ^ (4ri) -1 <5. Taking (11.21) and (11.51) 
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into account, in the case of integer l we have 

i 

|| x (>S' £ 6(D)Go) (•,i)||^« ( ^ ) ^ C'f ) ^c5- (I - s) ||(>S' £ 6(D)ao)(-,t)||^(o») 

s=0 
l 

^ci 6) ^5-^\\b(n)u 0 (-,t)\\ H ^on 

s =0 

l 

^ (da i) 1 ^ 2 ^<5 (/ 6) ||Duo(-,t)||^ S (o///). 

s=0 

(9.34) 


The norms of Duo(-,£) in 1*2(0) and /7 1 (0) have been estimated in (19.241) 
and (19.251) . By (|3.37l) (with O' replaced by O'"), we have 

l|Du 0 (-, t)\\ H s( 0 m^ ^ C' +1 4 s f 1 / 2 (8 2 + t l ) s / 2 e c ** /2 ||<p||l 2 (c»), s > 2. 

(9.35) 

From (19.241) . (19.251) . (19.341) . and (19.351) it follows that 

||x(5 £ 6(D)uo)(-,f)||^ (Rd) ^ C {8 h,- 1/2 (5~ 2 + t- 1 ) d/4 e- c * t/ ' 2 \\ip\\ L2{0) . 

(9.36) 

The constant C ® depends on d, «o, or, Ugllz,^, ||fl ,_1 ||Loo) an d the domain 
0. The case of half-integer l is considered similarly. 

As a result, relations (19.261) (19.291) . (19.331) . and (|9.36l) yield 

||/C D (f; e)(f - JC° D (t- e)(p\\ H i(o') < C (1) C (3) £ _1 e~ c * t/2 |M|l 2 ( 0 ) 

+ C'( 2 )(C'( 7 ) + C'(8))t- 1 /2( 5 -2 + 

This implies (13.381) with the constant C" = maxIC^C®; (C^ +0^)} 

and completes the proof of Lemma 3.15. 


9.9. Proof of Theorem 3.16. Inequality ((3.391) follows directly from 
(|3.36l) and (13.381) : herewith, = max{Ci6; Cn} + C^. 

Let us check ()3.40l) . From (13.391) and (1 1. 2[ ) . () 1,5 D it follows that 


\9 £ b(D) (e~ A ^ - e~ A ^ - e A e &(D) e -^) || L2(o) ^ 2(o0 

^ IblUoo (dai) 1/2 C d e/i d ((I; t)e~ c * t/2 . 
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We apply the identity f)9.171) . The norm of the second term in the right-hand 
side of (19.171) is estimated with the help of (1.5) and Lemma 9.1(1°): 

d 


e 


^ g £ b i A £ b j D i D j e ^ 
*,1=1 


L 2 (0)^L 2 (0') 


d 

< ll A£6 iX- D * L, i e "' 4 ^ii2(0)^i2(R d ) 

*,1=1 

d 

^ e llfflUoo Q; i(C ,(0 ^) 1 ^ 2 X] II xDiDje ADt \\L 2 (0)->H l - 1 (m.d-)-> 
*,1=1 


l = d/2. 

(9.39) 


Next, we apply Lemma 9.4. For the case of integer l we obtain 
d 

llxA D j e ADt \\L 2 {0)^H l ~ 1 {w i ) 

*.l=i 

^ dc\% yy 6 (l 1 v ’||e " 4 d4 ||l 2 (c>)-^//s+ 2 (c)//). 

s =o 


(9.40) 


The norm ||e ||l 2 (C))->h 2 (0) has been estimated in (13.51) . For s ^ 1, by 

(13.371) (with O' replaced by O"), we have 

l|e^||L 2 (0)^+ 2( 0") < C' s+ 2 2 s + 1 t“ 1 / 2 ( < r 2 + t-i)(^+i)/2 e -c*t/2_ ( 9 . 41 ) 

From (I9.39l) - (l9.4ip and (13.51) it follows that 


V 9 bi.\ bjDjl) 


—A 0 t 

e 


*,l=i 


L 2 (0)^L 2 (0') 


^ c ,(9 )er 1 / 2 (r 2 + r 1 )‘ i / 4 e- c * t/2 , 


(9.42) 


where the constant depends only on d, ao, ai, ||5 _1 ||l 00 , and the 

domain O. In the case of half-integer l inequality (19.42[) is checked similarly 
by using Lemma 9.4(2°). 

As a result, relations (19.171) . (|9.38l) . (|9.42[) . and (1.9) imply (13.401) with 
the constant = IbllLoo(dau) 1 / 2 ^ + C^\ This completes the proof of 
Theorem 3.16. 


9.10. Proof of Lemma 7.14. Arguments are similar to the proof of 
Lemma 3.15. As in Subsection 9.4, let p € 1^2(0; C n ) and let p = Vp. 
We put uo(-,f) = e - " 4 ^ 4 ^-), uo(-,f) = Pe>uo(-,t). According to (17.121) and 
(17.231) . IC]sr(t;£)ip = A e S e 6(D)uo and K.^(t;e)p = A e 6(D)uo- 
Similarly to (19.261) we obtain 


\\K. N (t-£)p - JC° N {t-£)p\\ H i io/) ^ ||A e x((5' e -/)6(D)u 0 )(-,t)||^i (Rd) . 
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Further considerations are similar to (19.271) - ()9. 371) : one should use estimates 
(17.41) . (17.51) . and (17.301) . This implies (17.311) . 

9.11. Proof of Theorem 7.15. Arguments are similar to the proof of 
Theorem 3.16. Inequality (17.321) follows directly from (I7.29|) and (17.311) . 
Herewith, % = max{Ci 6 ;Ci 7 } + 

Inequality (17.331) is deduced from (17.321) by analogy with (19.381) (19.421) . 
In addition, one should take (19.191) into account and apply estimates (17.51) 
and (17.301) . 
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